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Chapter  1 
Introduction 


The  impact  of  photonics  on  digital  communication  systems  is  extensive  and  well 
known.  Fiber  optics  carry  massive  amounts  of  data  between  users  and  services  around 
the  globe.  These  systems  are  finding  applications  in  shorter  and  shorter  distances, 
from  long-distance  telecommunications,  to  communication  between  servers  in  data 
centers,  to  interconnects  within  computers  themselves.  The  large  bandwidths  of  pho¬ 
tonic  systems  are  enabling  this  revolution. 

Less  well  known  are  the  benefits  of  photonics  to  high-frequency  analog  systems. 
These  “microwave  photonic”  systems  are  analogous  to  radio  systems,  where  baseband 
signals  are  modulated  onto  a  carrier  frequency.  Photonics  provide  very  high  frequency 
carriers,  around  194  THz  for  1550  nm  wavelength  light  used  with  standard  single  mode 
fiber,  so  signals  being  transmitted  and  manipulated  are  relatively  low  frequency  in 
comparison.  RADAR  and  wireless  communications  are  two  areas  that  can  greatly 
benefit  from  microwave  photonics. 

Improvements  in  the  analog  performance  of  photonic  systems,  especially  reduc¬ 
tions  in  noise  and  distortion,  have  direct  application  back  to  digital  communications. 
Next  generation,  commercial,  digital  fiber-optic  communication  systems  are  improv¬ 
ing  spectral  efficiency  (bits/s/Hz)  over  existing  fibers  in  order  to  save  on  infrastruc¬ 
ture  upgrades  to  fiber  optic  networks.  They  are  moving  away  from  simple  on-off-key 
(OOK)  representations  of  digital  data  in  favor  of  multi-level  and  coherent  modula¬ 
tion  techniques.  Optimizing  the  analog  performance  of  photonic  devices  and  systems 
increases  the  achievable  spectral  efficiencies  in  these  coherent  systems,  and  empowers 
this  next  advance  in  communications. 


1.1  Microwave  photonics  applications 

Microwave  photonics  is  the  study  of  photonic  devices,  such  as  lasers  and  photode¬ 
tectors,  performing  operations  at  microwave  frequencies,  and  the  application  of  these 
devices  to  microwave  systems.  Microwave  photonics  has  been  extensively  reviewed  by 
[1-4],  and  tutorial  information  has  been  published  in  book  form  by  [5,  6].  The  field 
broadly  defines  the  word  “microwave”  to  include  frequencies  ranging  from  hundreds 
of  megahertz  to  a  terahertz.  Much  work  has  been  performed  in  the  Super  High  Fre- 
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Figure  1.1:  Microwave  photonics  frequencies  of  interest. 
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quency  (SHF)  band,  defined  by  the  International  Telecommunication  Union  (ITU), 
which  ranges  from  3  GHz  to  30  GHz.  A  variety  of  RADAR  and  wireless  commu¬ 
nication  frequencies  fall  within  this  band.  Microwave  photonic  systems  are  analog 
systems.  They  are  analog  in  the  sense  that  they  manipulate  arbitrary  baseband  sig¬ 
nals  as  well  as  digital  signals  that  are  modulated  onto  a  higher  carrier  frequency. 

The  main  applications  for  microwave  photonics  can  be  categorized  into  signal 
transmission  and  signal  processing.  Photonics  can  be  used  for  antenna  remoting  and 
signal  distribution  for  a  variety  of  radio  technologies.  For  example,  an  array  of  CDMA 
antennas  are  used  to  extend  cellular  coverage  to  the  interior  of  a  large  building  such  as 
a  railway  station,  airport  or  subway.  Each  individual  antenna  transmits  the  detected 
signals  via  microwave  photonic  links  back  to  a  single  central  location  for  processing. 
With  the  right  design,  the  power  consumption  at  each  of  the  nodes  can  be  made  very 
small,  and  each  node  can  be  small  and  inexpensive  [2]. 

In  another  example,  [7],  an  array  of  radar  antennas  on  a  large  military  aircraft 
are  connected  to  a  central  location  with  microwave  photonic  links.  The  array  concept 
improves  the  overall  sensitivity  of  the  system  over  discrete  transmitters,  and  photonics 
allows  low-loss  collection  of  the  signals. 

Signal  processing  can  also  be  performed  with  microwave  photonics.  Researchers 
have  implemented  diverse  functions  such  as  tunable  bandpass  and  notch  filtering  of 
interference  [8],  microwave  mixing  [9],  arbitrary  waveform  generation  [10],  and  wide 
band  analog  to  digital  conversion  [11],  Photonics  can  be  used  for  the  generation 
of  microwave  signals.  Optoelectronic  oscillators  (OEOs)  are  one  technique  which 
can  produce  very  low-noise  microwave  oscillation  [12].  Photonics  can  also  generate 
millimeter  wave  signals  through  frequency  multiplication  techniques,  such  as  with 
with  injection  locked  lasers.  The  wide  bandwidth  of  microwave  photonics  makes  it 
ideal  for  performing  these  signal  processing  functions. 
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1.2  Advantages  for  signal  distribution 


For  signal  distribution,  the  competition  to  photonics  is  coaxial  cabling.  Conventional 
systems  are  fed  electronically  with  coaxial  cable  from  the  processing  station.  Elec¬ 
tronic  feeds  (which  are  3-300  meters  in  shipboard  and  avionics)  have  low  efficiencies  in 
size,  weight  and  power  (SWAP).  These  feeds  are  relatively  large,  inflexible  and  heavy 
because  of  multiple  coax  cable  runs.  They  have  high  loss,  which  limits  the  range  and 
requires  amplification  at  the  antennas.  Coax  is  not  especially  wide  bandwidth  because 
its  attenuation  is  frequency  dependent.  Coax  is  also  susceptible  to  electromagnetic 
interference  (EMI),  which  is  undesirable  in  military  applications. 

Microwave  photonic  links  have  been  explored  for  replacing  traditional  coaxial  links 
in  a  variety  of  applications  because  of  their  many  advantages  [13-15].  Optical  fibers 
have  significant  advantages  in  size  and  weight  over  microwave  coax.  Fiber  has  a  thin 
cross  section  and  its  bend  radius  is  much  tigher  than  for  coax.  By  remoting  signals 
with  fiber,  the  power  burden  can  be  shifted  to  the  processing  station.  Fibers  are  low 
loss,  and  the  loss  does  not  depend  very  much  on  the  signal  frequency.  Several  signals 
can  be  multiplexed  on  the  same  fiber  using  wavelength  division  multiplexing.  Fiber 
is  immune  from  EMI. 

The  most  successful  commercial  applications  have  been  in  hybrid-fiber-coax  (HFC) 
infrastructure  for  distributing  cable-television  signals  and  in  hybrid-fiber-radio  (HFR) 
for  distributing  cellular  signals  to  remote  antennas  [6,  13].  Military  radar  and  com¬ 
munication  systems  use  analog  fiber  optic  systems  for  antenna  remoting.  However, 
advanced  military  and  next  generation  wireless  systems  need  a  large  dynamic  range 
of  operation.  This  is  challenging  for  microwave  photonics,  as  they  are  not  yet  com¬ 
petitive  with  electronic  systems  in  terms  of  noise  and  distortion  [16].  In  addition, 
large  dynamic  range  is  important  for  microwave  photonics  signal  processing,  and  mi¬ 
crowave  photonic  links  are  a  performance  limiting  component  of  these  systems.  By 
improving  the  performance  of  the  microwave  photonic  links,  the  full  systems  also  are 
improved.  The  research  question  addressed  in  this  work  is  whether  we  can  have  the 
advantages  of  fiber  for  microwave  signal  transmission  while  still  maintaining  a  large 
dynamic  range. 

1.3  Dynamic  range  challenges 

The  dynamic  range  is  the  range  of  signal  amplitudes  that  can  be  transmitted  or  pro¬ 
cessed  by  a  system.  In  the  wireless  antenna  remoting  example,  the  dynamic  range  will 
play  a  role  in  determining  the  size  and  capacity  of  each  cell.  Remote-units  located 
close  to  the  antenna  have  to  limit  their  power  and  transmission  rate  if  they  exceed  the 
upper  end  of  the  dynamic  range,  and  remote  units  located  far  from  the  antenna  will 
not  be  noticeable  if  they  fall  below  the  lower  end  of  the  dynamic  range.  At  the  lower 
end,  the  range  is  limited  by  noise,  and  at  the  upper  end,  often  limited  by  the  point 
where  distortion  of  the  signal  by  the  system  is  noticeable.  Distortion  produces  har¬ 
monics  and  mixtures  between  signal  frequencies,  and  at  a  high  enough  signal  power, 
these  products  become  larger  than  the  noise.  This  particular  definition  of  dynamic 
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Figure  1.2:  Noise  and  distortion  limitations  on  the  dynamic  range  of  a  signal  trans¬ 
mission  system. 
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range  is  called  the  spurious  free  dynamic  range.  The  largest  distortion  products  tend 
to  be  the  second  and  third  orders,  which  grow  quadratically  and  cubically  with  the 
input  power.  Fig.  1.2  illustrates  the  concept  of  dynamic  range. 

In  the  following  sections,  I  will  define  relevant  concepts,  and  then  will  give  an 
example  comparison  between  an  electrical  link  and  a  microwave  photonic  link,  which 
shows  the  limitations  of  the  photonic  system  in  terms  of  dynamic  range. 

1.3.1  Microwave  photonic  links 

A  microwave  photonic  link  modulates  arbitrary  analog  signals  on  a  high  frequency 
carrier.  For  1550  nm  light,  the  carrier  is  approximately  194  THz.  The  analog  signals 
can  be  divided  into  frequency  bands,  for  example,  0.1-4  GHz,  4-8  GHz,  and  8-12 
GHz.  In  each,  an  RF  carrier  has  baseband  data  modulated  upon  it.  The  modulation 
process  creates  optical  sidebands  on  the  optical  carrier.  It  also  adds  noise  due  to  the 
phase  and  intensity  noise  of  the  laser,  and  distorts  the  signal.  The  detection  process 
recovers  the  electrical  signal,  but  also  adds  additional  noise  and  distortion  due  to 
shot  noise  and  nonlinearities  in  the  photodetection.  Fig.  1.4  illustrates  the  steps  in  a 
microwave  photonic  link. 

Typical  microwave  photonic  links  uses  intensity  modulation  and  direct  detection 
(IM-DD).  These  links  will  be  the  baseline  for  comparison  in  later  sections.  Fig.  1.4 
illustrates  a  direct  modulated  IM-DD  link,  where  the  bias  current  to  the  laser  is  varied 
with  the  signal,  thus  varying  the  intensity  of  the  emitted  light.  Fig.  1.5  illustrates  an 
externally  modulated  IM-DD  link,  where  a  lithium-niobate  Mach-Zehnder  modulator 
is  used  to  attenuate  the  laser  light  in  proportion  to  the  signal. 

1.3.2  Distortion 

Distortion  includes  both  harmonic  distortion  and  intermodulation  distortion.  Har¬ 
monic  distortion  creates  multiples  of  a  modulation  frequency.  It  is  typically  out-of- 
band,  but  this  is  still  important  for  multiband  links  and  ultra- wideband  links.  Inter¬ 
modulation  distortion  (IMD)  or  “intermod”  is  when  signals  of  different  frequencies  are 
mixed.  Typically,  the  most  important  IMD  terms  are  3rd  order  sum-and-difference 
products,  which  fall  in-band.  For  example,  for  two  modulation  frequencies  fi  and  /2, 
the  important  mixing  terms  are  2/2  —  fi  and  2/i  —  /2. 

Distortion  is  typically  quantified  using  a  two-tone-test.  Two  closely  spaced  fre- 
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Figure  1.3:  Diagram  of  signal  propagation  in  a  microwave  photonic  link.  The  output 
of  the  link  is  the  original  input  signal  with  the  addition  of  noise  and  distortion. 
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Figure  1.4:  Directly  modulated  IM-DD  link  comprised  of  a  semiconductor  laser,  op¬ 
tical  fiber  span  and  photodetector. 
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Figure  1.5:  Externally  modulated  IM-DD  link  comprised  of  a  laser,  Mach-Zehnder 

intensity  modulator,  optical  fiber  span,  and  photodetector. 
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Figure  1.6:  Harmonic  distortion. 
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Figure  1.7:  Intermodulation 
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quencies  are  transmited,  and  the  power  in  the  resulting  distortion  terms  are  measured 
with  a  spectrum  analyzer.  Interpolating  small  signal  measurements  to  high  input  pow¬ 
ers,  the  points  where  the  distortion  terms  are  equal  to  the  fundamental  in  power  are 
called  the  intercept  points.  The  output  powers  where  the  second-order  distortion  and 
third-order  distortion  are  expected  to  be  equal  to  the  fundamental  are  the  second- 
order  output  intercept  point  (OIP2)  and  third-order  output  intercept  point  (OIP3). 
Larger  values  for  OIP2  and  OIP3  mean  less  distortion. 

1.3.3  Noise 

Laser  relative  intensity  noise  (RIN),  laser  frequency  and  phase  noise,  optical  shot 
noise  and  modulator /detector  thermal  noise  all  contribute  to  the  noise  of  the  link. 
The  noise  of  the  link  is  quantified  by  its  noise  figure.  The  noise  figure  is  given  by 
the  input’s  signal  to  noise  ratio  divided  by  the  output’s  signal  to  noise  ratio,  usually 
assuming  the  input  is  thermal  noise  limited  in  a  50  ohm  impedance.  A  smaller  noise 
figure  link  introduces  less  noise.  The  noise  of  the  link  combined  with  the  distortion 
is  also  quantified  by  the  spurious  free  dynamic  range  (SFDR)  of  the  link. 
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Figure  1.8:  Two  tone  test. 
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Figure  1.9:  Output  intercept  points  and  spurious  free  dynamic  range. 


1.3.4  System  example 

I  would  like  to  give  an  example  that  illustrates  the  dynamic  range  of  a  very  good 
electronic  link  compared  to  a  microwave  photonic  link.  Suppose  I  have  to  transmit  a 
signal  centered  at  2  GHz  frequency  over  a  distance  of  100  m.  Very  low  attenuation, 
high  performance  coaxial  cabling  has  been  developed  for  avionics.  At  best,  these 
cables  have  an  attenuation  of  20  dB  per  100  m.  Typical  commercial  cabling  has  much 
higher  attenuation. 

Assume  I  place  a  high-dynamic-range  pre-amplifier  before  the  link  to  overcome 
the  20  dB  attenuation.  I  assume  a  gain  of  20  dB,  a  1  dB  noise  figure,  and  a  third- 
order  output  intercept  point  of  10  W  (40  dBm).  Amplifiers  are  typically  limited  by 
third-order  distortion,  so  the  OIP3  value  is  relevant  to  calculating  the  spurious  free 
dynamic  range.  In  decibel  units,  the  SFDR  is  given  by 


SFDR 


\  ( OIP3  -  G  +  174  ^  -  10 log10  (B)  -  NF)  . 


where  G  is  the  gain  in  dB  units  and  B  is  the  bandwidth.  In  1  Hz  bandwidth,  this 
would  give  a  dynamic  range  of  129  dB.  (75  dB  in  100  MHz  of  bandwidth).  The  link 
noise  figure  is  limited  to  the  noise  figure  of  the  amplifier,  and  is  about  1  dB. 

I  will  next  illustrate  the  dynamic-range  of  a  typical  photonic  link  using  commer¬ 
cially  available  components.  This  system  consists  of  an  electrical  to  optical  (e-to-o) 
transducer,  a  fiber  optic  transmission  line,  and  an  optical  to  electrical  (o-to-e)  trans- 
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Figure  1.10:  Electrical  link  dynamic  range  example  with  list  of  typical  parameters. 


Parameter 

Value 

Signal  frequency 

2  GHz 

Distance 

100  m 

Coaxial  cable 

Low  loss  PTFE  dielectric  or  0.325  in  rigid  coax 

Attenuation 

20  dB  /  100  m 

Amplifier  gain 

20  dB 

Amplifier  noise  figure 

1  clB 

Amplifier  OIP3 

10  W  (40  dBm) 

Spurious  free  dynamic  range 

129  dB  in  1  Hz  bandwidth 

Noise  figure 

1  dB 

ducer.  Our  e-to-o  transducer  is  a  high  efficiency  Mach-Zehnder  modulator,  which 
modulates  a  microwave  signal  onto  the  intensity  of  an  optical  carrier  provided  by 
a  semiconductor  laser.  The  o-to-e  transducer  is  a  photodiode,  which  detects  the 
envelope  of  the  intensity  modulation.  For  100  m  of  single- mode  optical  fiber,  the 
transmission  loss  is  less  than  0.05  dB,  which  is  why  fiber  optics  are  extensively  used 
for  long  distance  communications.  The  parameters  below  were  chosen  to  give  a  gain 
of  0  dB  for  the  link. 

The  e-to-o  transducer  has  a  sinusoidal  transfer  function  of  light  intensity  versus 
voltage,  which  contributes  a  large  amount  of  distortion  to  the  final  signal.  This  system 
requires  a  photodiode  capable  of  handling  high  optical  power.  Research  devices  have 
been  demonstrated  that  can  handle  much  higher  powers  than  this,  but  this  is  still 
an  expensive  device.  The  third-order  distortion  and  shot  noise  limited  SFDR  for  this 
link  is  derived  in  dB  units  per  1  Hz  bandwidth  by  [17]  as 


SFDR 


(1.1) 


where  e  is  the  elementary  charge  and  Idc  the  effective  DC  photocurrent.  In  1  Hz  of 
bandwidth,  this  would  give  a  dynamic  range  of  116  dB,  which  is  13  dB  worse  than 
the  electronics  case.  What’s  worse  is  the  noise  figure  of  this  particular  link,  which  is 
18.5  dB,  compared  to  1  dB  for  the  electronics  case.  Assuming  a  shot-noise  limited 
receiver,  the  noise  figure  is  calculated  by  using  [17] 

2eK2 


NF  =  10  log10 


Id.cF2KT  Zin 


(1.2) 


Figure  1.11:  Photonic  link  dynamic  range  example  with  list  of  typical  parameters. 
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Parameter 

Value 

Signal  frequency 

2  GHz 

Distance 

100  m 

Fiber  attenuation 

<  0.05  dB 

Modulator 

High  efficiency  Lithium  Niobate  MZM 

Halfwave  voltage 

3  V 

Photodetector 

High  power  InGaAs  PIN  photodiode 

Photocurrent 

20  mA 

Spurious  free  dynamic  range 

116  dB  in  1  Hz  bandwidth 

Noise  figure 

18.5  dB 

where  Vn  is  the  modulator  half-wave  voltage,  K  is  Boltzmann’s  constant,  T  is  the 
system  temperature  (300  K),  Zin  is  the  input  impedance  of  the  system,  typically  50 
ohms. 

The  noise  figure  is  heavily  influenced  by  the  inefficiency  of  the  e-to-o  transducer, 
given  by  large  Vn  .  In  addition,  in  a  real  system,  the  input  and  output  of  the  system 
must  be  impedance  matched.  If  passive  impedance  matching  is  used,  the  usable  signal 
level  is  further  reduced.  For  better  noise  and  dynamic  range  performance,  I  would 
like  to  have  higher  efficiency  e-to-o  conversion,  and  e-to-o  conversion  that  is  much 
more  linear. 


1.4  Techniques  to  improve  dynamic  range 

There  has  been  much  work  performed  to  improve  the  dynamic  range  of  microwave 
photonic  links  through  both  optical  design  and  by  using  electrical  system  techniques. 
The  noise  and  linearity  performance  of  externally  modulated  photonic  links  scale 
with  increasing  optical  power  at  the  detector,  as  can  be  seen  in  equations  1.1  and  1.2. 
Work  has  been  dedicated  to  improving  the  power  handling  of  photodetectors  and 
their  linearity  [18,  19],  designing  high  power  handling  optical  fibers  to  reduce  optical 
power  induced  stimulated  Brillouin  scattering,  and  reducing  laser  relative  intensity 
noise  to  ensure  that  the  receiver  is  shot  noise  limited  at  higher  optical  powers.  On  the 
modulator  side,  there  have  been  efforts  to  decrease  the  halfwave  voltages  of  Mach- 
Zehnder  modulators  to  improve  the  link  gain. 
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Researchers  have  developed  modulator  designs  which  improve  the  link  linearity 
over  that  of  a  simple  MZM.  These  modulators,  with  multiple  modulation  sections, 
have  a  transfer  function  that  is  more  linear  than  the  MZM’s  sinusoidal  one  [20]. 
However,  linearized  modulators  are  complicated,  difficult  to  fabricate,  difficult  to  op¬ 
timize  for  high-frequency  (traveling-wave)  operation,  and  have  had  little  experimental 
demonstration. 

Laser  designers  have  worked  on  improving  the  direct  intensity  modulation  linearity 
of  semiconductor  lasers.  There  has  been  interest  in  modeling  and  choosing  physical 
device  parameters  which  minimize  the  distortion  (for  example,  [21]).  Strong  optical 
injection  locking  is  one  technique  which  has  been  shown  improve  to  linearity  by  in¬ 
creasing  laser  resonance  frequency  [22],  System  design  techniques,  including  using  a 
push-pull  configuration  with  balanced  detection  have  shown  some  success  [15]. 

There  are  electronic  means  for  improving  link  distortion  by  compensating  for  mod¬ 
ulation  nonlinearity.  These  include  predistortion  [23,  24],  feedforward  linearization 
techniques  [25],  and  feedback  linearization  [26].  However,  these  techniques  require  fast 
electronics  to  perform  the  linearization.  At  the  present  time,  they  are  not  useable  for 
very  high  frequency  microwave  photonics  beyond  a  few  GHz. 

In  this  work,  I  have  demonstrated  linearity  improvement  using  two  techniques 
called  phase  modulation  direct-detection  (PM-DD)  and  frequency  modulation  direct- 
detection  (FM-DD).  These  approaches  are  based  on  optical  system  design  and  do 
not  require  high-speed  electronics  for  linearization,  so  they  are  potentially  useable  to 
very  high  modulation  frequencies.  The  modulation  techniques  are  simple,  requiring 
only  a  lithium  niobate  phase  modulator  or  a  direct-modulated  multi-section  laser. 
The  demodulation  process  does  require  optical  filters,  but  these  are  realizeable  with  a 
variety  of  fabrication  technologies.  PM-DD  and  FM-DD  systems  scale  in  performance 
with  detector  power  handling  as  do  IM-DD  links,  so  they  benefit  from  more  general 
device  research  in  the  field.  The  following  chapters  will  present  theoretical  derivations, 
simulations  and  experimental  evidence  of  the  benefits  which  PM-DD  and  FM-DD 
microwave  photonic  links  can  provide  to  improve  the  noise  and  linearity  in  microwave 
photonic  systems. 
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Chapter  2 


Theory  of  PM-DD  and  FM-DD 
links 

2.1  Motivation  for  phase  and  frequency 
modulation 

Microwave  photonic  links  (MPLs)  with  large  dynamic  range  are  an  essential  com¬ 
ponent  of  high-performance  microwave  distribution  and  processing  systems.  Large 
dynamic  ranges  require  low  signal  distortion  and  low  noise  figures.  These  metrics  are 
poor  in  traditional  intensity  modulated  links,  but  modulation  is  not  limited  to  the 
intensity.  Other  parameters  of  the  light  can  be  used  to  convey  information,  including 
the  amplitude,  phase,  frequency,  spatial  modes,  and  polarization  of  the  light’s  electric 
held.  Phase  modulation  (PM)  and  frequency  modulation  (FM),  where  the  instaneous 
optical  phase  or  frequency  is  varied  in  proportion  to  the  input  signal,  are  consid¬ 
ered  to  be  promising  alternatives  to  IM.  PM  is  a  promising  modulation  technique 
for  MPLs  because  devices  are  highly  linear.  Phase  modulators  based  on  the  linear 
electro-optic  effect,  including  those  fabricated  in  lithium  niobate,  are  intrinsically  lin¬ 
ear,  and  authors  have  also  reported  linear,  integrable  phase  modulators  fabricated  in 
indium-phosphide  [27]. 

The  signal  loss  of  MPLs  is  an  important  factor  for  links  and  systems  as  it  impacts 
the  signal  to  noise  ratio.  Traditional  intensity-modulated  direct-detection  (IMDD) 
links  experience  large  signal-loss  and  resulting  low  noise  figures  due  to  the  low  modu¬ 
lation  efficiency  of  lithium  niobate  Mach  Zehnder  modulators  (MZMs).  On  the  other 
hand,  directly  modulated  frequency  modulated  (FM)  lasers  have  been  demonstrated 
with  high  modulation  efficiency  and  with  modulation  bandwidths  that  are  not  limited 
by  the  laser  relaxation  frequency  [28].  Recent  work  on  multi-section  DFB  [29]  and 
EML  lasers  [30]  have  produced  modulation  efficiencies  two  orders  of  magnitude  better 
than  traditional  intensity  modulation.  An  improvement  in  modulation  efficiency  could 
make  a  major  impact  on  the  noise  performance  of  microwave  photonic  links.  Besides 
high  modulation  efficiency,  the  performance  of  these  devices  is  also  more  linear  than 
direct  intensity  modulation  and  Mach  Zehnder  modulators,  and  there  is  low  thermal 
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Figure  2.1:  Externally  modulated  PM-DD  link  comprised  of  a  laser,  lithium  niobate 

phase  modulator,  optical  fiber  span,  optical  filters  and  photodetector. 
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Figure  2.2:  Directly  modulated  FM-DD  link  comprised  of  a  multi-section  laser,  optical 

fiber  span,  optical  filters  and  photodetector. 

Gain  Bias  RF  In 


cross-talk  in  integrated  laser  arrays.  PM  and  FM  have  favorable  characteristics  for 
linearity  and  gain  in  MPLs. 

2.2  Link  architecture 

Because  photodiodes  only  respond  to  the  intensity  envelope  of  the  light,  phase  and 
frequency  modulation  can  not  be  directly  detected.  Coherent  detection  using  hetero¬ 
dyning  is  one  possibilc  demodulation  scheme,  but  heterodyning  is  nonlinear  and  adds 
complexity.  Alternatively,  one  can  use  a  direct-detection  system.  We  have  designed 
demodulators  which  use  optical  filters  to  convert  the  phase  and  frequency  modulation 
into  AM  before  direct  detection  at  a  photodetector.  The  filters  are  called  phase  and 
frequency  discriminators.  The  demodulation  process  is  called  phase-modulation  or 
frequency-modulation  direct-detection  (PM-DD  or  FM-DD  [31]),  filter-slope  detec¬ 
tion,  or  interferometric  detection  [17].  The  architecture  for  the  PM-DD  and  FM-DD 
links  consists  of  a  modulation  source,  discriminator  filters,  and  balanced  detectors. 
The  link  architectures  are  shown  in  Fig.  2.1  and  Fig.  2.2.  Discriminators  for  PM- 
DD  and  FM-DD  links  have  similar  design  because  PM  is  identical  to  FM  but  with  a 
modulation  depth  that  is  linearly  dependent  on  modulation-frequency. 

The  sidebands  of  a  phase-modulated  or  frequency-modulated  signal  possess  certain 
amplitude  and  phase  relationships  among  themselves  such  that  the  envelope  of  the 
signal  is  independent  of  time.  A  discriminator  works  by  modifying  these  phase  and 
amplitude  relationships  such  that  the  amplitude  of  the  envelope  of  the  resultant  signal 
fluctuates  in  the  same  manner  versus  time  as  did  the  instantaneous  frequency  of  the 
original  signal  [32],  One  can  also  think  of  the  discriminator  as  a  filtering  function  with 
a  frequency  dependent  amplitude.  The  slope  of  the  function  converts  variations  in 
the  optical  frequency  into  variations  in  the  amplitude.  This  view  is  accurate  for  slow 
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Figure  2.3:  PM-DD  link  using  a  Mach  Zehnder  interferometer,  and  an  IM-DD  link 
with  a  dual-output  Mach  Zehnder  modulator.  For  a  given  photocurrent,  these  links 
have  the  same  figures  of  merit.  The  IM-DD  link  may  use  a  multiplexing  scheme  to 

combine  both  complementary  signals  onto  the  same  optical  fiber. 
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RF  Out 
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variations  of  the  optical  frequency.  However,  it  generally  can  be  misleading  since  it 
assumes  that  the  instantaneous  frequency  of  the  light  is  equivalent  to  a  time-averaged 
frequency.  Nevertheless,  the  model  is  instructive  as  it  suggests  that  functions  with 
larger  slopes  will  have  higher  conversion  efficiency  to  AM,  and  that  a  function  with 
many  large  high  order  derivatives  will  distort  the  AM  signal  more  than  one  with  a 
more  “linear”  function. 

The  system’s  performance  is  determined  by  the  transfer  function  of  the  optical 
filter.  For  example,  a  Mach  Zehnder  interferometer  (MZI)  after  a  phase  modulator 
has  comparable  nonlinearity  to  a  Mach  Zehnder  modulator  [17].  This  is  shown  in  Fig. 
2.3.  Authors  have  proposed  various  discriminator-filters  to  optimize  the  demodulation 
for  low  distortion,  including  birefringent  crystals  [33],  asymmetrical  Mach  Zehnder 
interferometers  (a-MZI)  [17,  34],  Fabry-Perot  filters  [35],  fiber  Bragg  gratings  [36] 
and  tunable  integrated  filters  [37,  38]. 

In  the  PM-DD  and  FM-DD  links,  the  ideal  transfer  function  of  the  optical  filter 
is  a  linear  ramp  of  field-transmission  versus  offset  frequency  from  the  optical  carrier, 
which  is  a  quadratic  ramp  of  power  transmission.  The  ideal  filters  have  linear  phase. 
The  power  is  split  between  two  filters  with  complementary  slope,  and  detected  with  a 
balanced  photodetector.  I  first  analyzed  this  complementary  linear-field  demodulation 
scheme  analytically  in  [39].  The  link  architecture  is  shown  in  Fig.  2.4.  A  single 
filter  and  detector  has  low  third-order  distortion,  and  the  balanced  detection  cancels 
second-harmonics  of  the  signal’s  Fourier-frequency  components  produced  by  squaring 
of  the  AM.  Since  it  is  difficult  to  implement  this  transfer  function  in  optics,  a  realized 
discriminator  will  have  a  transfer  function  with  some  non-idealities. 
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Figure  2.4:  Ideal  filter  transfer  functions  for  an  optical  PM  or  FM  discriminator  in  a 

complementary  linear-field  demodulation  scheme. 
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2.3  History 

The  work  of  Harris,  [40],  was  the  earliest  use  of  a  quadrature  biased  Mach  Zehnder 
interferometer  structure  to  discriminate  optical  FM.  An  interferometric  path  differ¬ 
ence  was  created  by  passing  the  light  through  a  birefringent  crystal  when  the  light’s 
polarization  was  angled  between  the  fast  and  slow  axes  of  the  crystal.  It  was  noted 
by  Harris  that  optimal  FM  to  AM  conversion  occurs  at  the  quadrature  bias  point. 
The  technique  was  also  applied  to  phase  modulated  light  in  [32] .  Besides  PM  to  AM 
discrimination,  suppression  of  unwanted  incident  AM  was  done  by  applying  a  180 
degree  phase  shift  to  one  of  the  two  complementary  polarization  states  at  the  output 
of  the  discriminator.  The  initial  AM  canceled  when  both  polarization  states,  now 
with  their  PM  in  phase  but  AM  180  degrees  out  of  phase,  were  detected  at  a  single 
polarization-insensitive  photodetector. 

Another  physical  implementation  of  the  MZI  style  discriminator  using  mirrors  and 
beam  splitters  was  suggested  by  [34],  In  this  case,  balanced  photodetection  was  used 
to  cancel  AM.  Such  an  interferometer  was  experimentally  verified  by  [41].  [34]  also 
suggested  the  use  of  balanced  detection  for  the  birefringent  crystal  device  of  [40] . 

Concurrent  to  the  development  of  direct  frequency  modulation  of  semiconductor 
lasers  in  works  such  as  [42],  [43]  performed  digital  data  transmission  experiments  using 
a  Michelson  interferometer  to  discriminate  optical  frequency  shift  keying  (FSK). 

The  use  of  FM  semiconductor  lasers  and  discriminator  detection  was  extended  to 
transmitting  subcarrier-multiplexed,  analog  signals  for  applications  in  cable  television 
distribution.  Experimental  results  for  a  Fabry-Perot  discriminated,  FM  subcarrier- 
multiplexed  system  were  presented  by  [44] .  An  array  of  optical  frequency  modulated 
DFB  lasers  and  a  Fabry-Perot  discriminator  were  used  to  transmit  and  demodulate 
a  large  number  of  microwave  FM,  analog  video  channels.  A  similar  system  was  also 
used  to  transmit  subcarrier-multiplexed,  digital  signals  in  [35]. 

Because  analog  links  require  high  linearity  and  low  noise,  a  number  of  authors  have 
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derived  figures  of  merit  for  the  performance  of  analog  FM-DD  links.  [45]  analyzed  the 
frequency-dependent  response  of  a  link  with  a  quadrature  biased  MZI  discriminator 
subject  to  large  modulation-depth  AM  and  FM.  [46]  studied  the  intermodulation 
distortion  for  a  Fabry-Perot  discriminated  link  with  a  large  number  of  channels,  while 
taking  into  account  both  FM  and  IM  on  each  channel.  [17]  derived  figures  of  merit  for 
the  dynamic  range  of  a  phase  modulated  link  with  an  MZI  discriminator  and  balanced 
detection. 

[47]  studied  a  link  with  an  arbitrary  discriminator.  The  general  formulae  were 
applied  to  the  particular  cases  of  an  MZI  and  a  Fabry-Perot  interferometer.  How¬ 
ever,  the  analysis  was  inaccurate  since  it  looked  at  the  system  in  terms  of  light  in¬ 
tensity  transmission  through  the  interferometer,  and  ignored  the  coherence  of  the 
filtering.  The  transmission  was  expanded  in  terms  of  a  Taylor  series.  The  analysis 
assumed  that  derivatives  of  the  transmission  spectrum  of  the  interferometer  (in  the 
Fourier-frequency  domain)  with  respect  to  the  instantaneous  optical  frequency  were 
proportional  to  overall  link  nonlinearity.  Similar  (inaccurate)  theoretical  analyses  us¬ 
ing  Taylor  series  were  published  by  [48]  and  [49].  However,  these  papers  did  include 
new  models  for  the  nonlinearities  in  the  lasers’  FM  and  included  the  effects  of  residual 
IM. 

To  improve  the  linearity  of  an  FM-DD  link,  many  alternatives  to  the  Mach- 
Zehnder  and  Fabry-Perot  interferometers  have  been  suggested.  In  very  early  work, 
[33]  proposed  a  linear-field  discriminator  using  a  network  of  birefringent  crystals.  The 
device  was  a  tenth-order  finite-impulse-response  (FIR)  filter.  The  series  of  crystals 
worked  as  a  series  of  cascaded  Mach  Zehnder  interferometers,  and  the  network  was 
equivalent  to  a  lattice  filter  architecture.  The  filter  coefficients  chosen  were  the  ex¬ 
ponential  Fourier  series  approximation  to  a  triangular  wave.  The  authors  understood 
that  linear  demodulation,  required  for  high  fidelity  signal  transmission,  could  be  ac¬ 
complished  with  a  discriminator  that  has  a  linear  FM  to  AM  transfer  function,  and 
that  high-order  filters  could  be  used  to  implement  this  linear-field  discriminator. 

Except  for  the  early  work  of  [33] ,  other  “linearized”  discriminators  in  the  literature 
were  designed  such  that  the  filter’s  optical  intensity  transmission  ramped  linearly  with 
frequency  offset  from  the  carrier,  rather  than  the  field  amplitude.  These  designs  are 
not  consistent  with  our  theoretical  link  models.  [50]  and  [51]  proposed  pairs  of  chirped 
fiber-Bragg  gratings  with  either  the  index  variation  or  chirp  rate  varied  nonlinearly. 
[38]  proposed  a  frequency  discriminator  based  on  an  MZI  with  ring  resonators  in  its 
arms.  [52]  suggested  that  the  linearity  of  a  Sagnac  discriminator  could  be  improved 
by  adding  ring  resonators. 

There  have  been  recent  experimental  results  for  discriminators  with  intensity  ver¬ 
sus  frequency  offset  ramps.  None  of  these  devices  have  demonstrated  significant 
linearity  improvement  over  a  MZI  .  Design  and  experimental  results  from  a  micro¬ 
ring  structure  implemented  in  a  CMOS  waveguide  process  were  reported  by  [37,  53]. 
Experimental  and  theoretical  results  using  fiber-Bragg  gratings  were  presented  in 
[36,  54-58].  These  experiments  used  pairs  of  complementary  gratings  designed  to 
have  a  a  transfer  function  whose  intensity  transmission  ramped  linearly  with  offset 
frequency  from  the  carrier.  The  gratings  were  low-biased  to  perform  carrier  suppres¬ 
sion.  In  [56,  58],  the  authors  presented  a  clipping- free  dynamic  range  limit  for  an 
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FM-DD  system.  (In  related  work,  [59,  60],  the  authors  used  Bragg  gratings  to  con¬ 
vert  phase  modulation  into  single  sideband  modulation.)  After  a  theoretical  analysis, 
the  authors  later  realized  the  limitations  of  their  discriminator  filter  design,  [57]: 

[...]  the  ideal  linear  power  reflectivity- versus- frequency  curve  does  not 
result  in  an  ideal  half-wave  rectification,  as  suggested  by  the  simple 
time-domain  view.  Rather,  in  addition  to  the  signal  component,  the 
output  includes  a  dc  component  as  well  as  a  nonlinear  distortion. 

They  explained  the  discrepancy,  [36]: 

The  reason  this  intuition  fails  is  that  combining  a  time-domain  view  of 
the  FM  signal  (instantaneous  frequency,  not  averaged  over  time)  with 
a  frequency  domain  view  of  the  FBG  filter  response  is  inconsistent 
with  the  frequency  domain  analysis  [...] 

It  is  erroneous  to  think  of  the  modulated  signal  in  terms  of  its  instantaneous  frequency 
while  looking  at  the  frequency  spectrum  of  the  filter.  The  carrier  is  not  really  being 
swept  along  the  ramp  of  the  filter  by  the  modulation,  so  considering  it  in  the  same 
way  as,  for  example,  the  small-signal  current  to  voltage  relationship  of  an  amplifier 
is  not  correct.  In  this  work,  I  present  complementary  linear-field  demodulation  as  a 
technique  that  can  produce  a  microwave  photonic  link  with  low  distortion. 

2.4  Analytical  link  analysis 

In  this  section,  I  derive  hgures-of-merit  for  a  PM-DD  or  an  FM-DD  link  that  uses 
an  arbitrary  optical  filter  for  discrimination,  following  my  published  work  in  [39]. 
This  derivation  is  related  to  earlier  theoretical  work  by  [36] ,  who  published  results  for 
single-tone  modulation.  Follow-up  work  has  been  performed  by  [61],  which  consider 
links  with  partially  coherent  sources.  I  find  expressions  for  the  currents  at  each 
microwave  frequency  at  the  output  of  the  link  under  a  two-tone  test.  I  take  a  small- 
modulation-depth  approximation.  The  standard  definitions  for  the  linearity  figures 
of  merit  rely  on  this  small  signal  approximation.  I  obtain  expressions  for  the  signal- 
to-noise  ration  (SNR),  second-order  and  third-order  output  intercept  points  (OIP2 
and  OIP3),  spurious-free  dynamic  range  (SFDR)  and  noise  figure  (NF).  I  apply  these 
general  formulae  to  the  specific  cases  of  the  Mach  Zehnder  interferometer,  a  linear 
intensity  ramp  filter  and  complementary  linear-field  filters.  For  the  linear-field  filter, 
I  derive  the  noise  figure’s  dependence  on  the  link’s  regime  of  operation  and  quantify 
the  effect  of  filter  curvature  and  the  laser’s  residual  IM  on  the  distortion. 

2.4.1  Two  tone  derivation 

An  optical  signal  that  is  phase  or  frequency  modulated  by  two  sinusoidal  tones  can 
be  represented  by  the  time  varying  electric  held 

emod  (t)  =  Ky/2Popt  cos  [27r fct  +  /3i  sin  (27 r/R)  +  /32  sin  (2n f2t)]  (2.1) 
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where  Popt  is  the  rms  optical  power,  k  is  a  constant  with  units  relating  optical  held 
and  optical  power  such  that  Popt  =  (e  (f)2)  /k2,  fc  is  the  frequency  of  the  optical 
carrier,  fi  and  f2  are  the  modulation  frequencies  and  Pi  and  p2  are  the  angle  mod¬ 
ulation  depths.  For  PM,  P  is  the  peak  phase  shift  induced  by  the  modulator.  For  a 
peak  applied  voltage  of  V,  the  peak  phase  shift  is  j3  —  7rV/Vn  (/),  and  the  halfway 
voltage  is  generally  frequency  dependent.  For  FM,  each  modulation  depth  is  equal  to 
the  maximum  optical  frequency  deviation  of  the  carrier  induced  by  the  modulation 
divided  by  the  frequency  of  the  modulation,  f3  =  Sf/f  .  The  modulation  of  the  light 
can  be  thought  of  in  terms  of  variations  in  the  instantaneous  frequency  of  the  light 
due  to  the  applied  signal.  The  optical  frequency,  or  wavelength,  varies  sinusoidally  in 
time.  The  instantaneous  frequency  of  the  light  is  given  by  the  derivative  of  the  phase 
of  the  light, 

77“  [2vr fct  +  Pi  sin  (2 7t/R)  +  P2  sin  (27r/2t)]  =  fc  +  Sfl  cos  (2 nfit)  +  5 h  cos  (2vr  f2t) 

2  n  at 

(2*2) 

The  link  generally  has  additional  undesired  residual  IM  and  noise.  The  correction 
to  the  electric  held  is 

emod  ( t )  —a  ( t )  +  K^2Popt  [1  +  n  (£)]  (2.3) 

•  \Jl  +  mi  cos  (2nfit  +  (f>)  +  m2  cos  JZnfet  +  0) 

•  cos  [2nfct  +  Pi  sin  (2nfit)  +  p2  sin  (2nf2t)  +  <p{t)\ 

where  n  (t)  is  the  RIN  of  the  source,  (p(t)  is  the  phase  noise  of  the  source,  a  (t)  is  the 
ASE  noise  from  an  optical  amplifier,  mi  and  m2  represent  the  IM  depths  for  the  two 
tones  and  p  is  the  phase  difference  between  the  IM  and  the  FM.  The  link  will  also 
amplify  thermal  noise  present  at  the  input. 

In  the  next  few  equations,  I  expand  the  expression  for  the  modulated  electric 
held  into  its  frequency  components  so  that  filtering  can  be  expressed  in  the  frequency 
domain.  The  residual  IM  depth  and  the  intensity  noise  are  assumed  to  be  much 
smaller  than  the  angle  modulation,  so  the  square  root  in  (2.3)  can  be  expanded  using 
a  Taylor  series,  yielding 

emod  (' t )  «a  (t)  +  K^/2Popt  (2.4) 

1  +  cos  (27r fit  +  p)  +  ^m2  cos  ( 2irf2t  +  p)  +  )ai  (t) 

■  cos  [27r/ct  +  Pi  sin  (2-nfit)  +  p2  sin  (2ir  f2t)  +  ip(t)\ 

Ignoring  noise,  this  can  be  written  using  an  angular  addition  trigonometric  identity 
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as 


^ mod  (t")  ^  \/2  r<>Pl  Av  ^ 

cos  [2nfct  +  Pi  sin  {2n fit)  +  P2  sin  (2n  f2t)] 

+  cos  [27 r  {fc  +  fi)t  +  Pi  sin  ( 2nfit )  +  p2  sin  {2nf2t)  +  0] 

+  cos  [2tt  (/c  -  /i)  t  +  Pi  sin  (2v ifp)  +  p2  sin  (2vr  f2t)  -  p\ 

+  Kn2  cos  [2vr  (/c  +  f2)  t  +  ft  sin  (2v ifp)  +  p2  sin  (2vr  f2t)  +  0] 

+  ^-m2  cos  [2vr  (/c  -  /i)  t  +  ft  sin  (2v t/^)  +  p2  sin  (: 2nf2t )  -  0]  | 

The  Jacobi-Anger  expansion  is  given  by  el^cose  =  '}^=_00jnJn  (P)  em6 where  j  is  the 
imaginary  unit  and  ft  (ft  is  a  Bessel  function  of  the  first  kind.  Applying  this  formula, 
the  expression  in  final  form  expands  to 


£ mod  (t)  1^  V2P '}Pt 

oo  oo 

55  55  ft  (ft)  ft  (ft)  exP  b’27r  (ft  +  ™/i  +pf2)t] 

n=— oo  p=—oo 

^  oo  oo 

+  -mi  E  E  ft  {Pi)  ft  (ft)  exp  [j2vr  (ft  +  [n  +  l]/i  +  pf2)  t  +  j(j)\ 

n=— 00  p=—oo 
00  00 

+  4mi  51  51  ft  (ft)  ft  (ft)  exp  [j2vr  (ft  +  [n  -  l]ft  +  pft)  f  -  jft 

n= — 00  p=—c*o 
^  00  00 

+  -ra2  E  E  ft  {Pi)  Jp  {P2)  exp  [j27T  (ft  +  nfi  +  [p  +  l]/2)  t  +  jft] 

n= — 00  p=—oo 
1  00  00 

+  4m2  55  55  (ft)  Jp  (ft)  exp  [j2vr  (ft  +  n/,  +  [p  -  l]/2)  t  -  .ft]  l 

n=—oo  p=— 00  J 

An  arbitrary  optical  filter  is  used  on  the  link  to  convert  the  angle  modulation 
to  IM.  With  multiple  detectors,  we  denote  the  transfer  function  seen  by  the  field 
before  each  detector  as  Hz  (/)  for  the  zth  of  Z  detectors.  For  example,  Fft  (ft)  is  the 
attenuation  of  the  optical  carrier  seen  at  the  first  detector.  The  transfer  function 
includes  the  splitting  loss.  For  later  convenience,  I  employ  a  shorthand  notation  for 
electric  field  transmission  at  each  frequency  in  the  optical  spectrum  that  corresponds 
to  an  optical  sideband  around  the  carrier: 


hn,P  =  Hz  (ft  +  nfi  +  pf2)  (2.5) 

where  n  and  p  are  integer  indices  and  H  is  the  complex  transfer  function  of  the 
filter,  representing  its  phase  and  amplitude  response,  including  any  insertion  losses 
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or  optical  amplifier  gain.  For  example,  /i0,o  is  the  field  transmission  for  the  optical 
carrier,  and  /r_i,o  is  the  transmission  of  the  negative,  first  order  sideband  spaced  ft 
away  from  the  carrier. 

The  electric  field  after  the  filter  at  photodetector  z  is 

edet  (t)  =emod  (t)  *  hZ  (t) 

^  \/^  F*opt 
oo  oo 

J”(ft)  ft  (ft)  K,p  exp  b'27r(/c  +  n/i  +pf2)t] 

n=— oo  p=— oo 

^  oo  oo 

+  -mi  E  E  Jn  (ft)  JP  (ft)  hzn+l  p  exp  \j2n  (ft  +  [n  +  l]/i  +  pf2)  t  +  j<j>\ 

n= — oo  p=—oo 
1  oo  oo 

+  4mi  5Z  5Z  Jn(/3i)4(/?2)^_i)Pexp[j27r(/c+[n-l]/i+p/2)i-i(/)] 

n=—oo  p=— oo 
1  X  x 

+  -m2  E  E  Jn  (ft)  (ft)  ^n,p+l  eXP  [j'27r  ( fc  +  nfi  +  \p+l]f2)t  +  j(f>] 


n=— 00  p=—oo 
00  00 


+  4m2  X]  Jn  (ft)  ft  (ft)  ^n.p-1  eXP  h'27r  (fc  +  nfi  +  \p-  1  ]/2)  t 

Tl— —  OO  p=  —  00 

The  indices  of  each  infinite  sum  can  be  renumbered  to  obtain 


edet  (t)  —  K^/2PoptRe< 

OO  OO 

E  E  Jn  (ft)  Jp  (ft)  KP  exp  [j2vr  (ft  +  nft  +  pf2)  t] 

:  —  O 

1 


n=—oo  p=— 00 

00  00 


mi  E  E  ft-i  (ft)  ft  (ft)  exp  [j'27t  (ft  +  nfi  +  pf2)  t  +  j<j>\ 


Tl— —  OO  p=— OO 

00  00 


mi  E  E  d re+ 1  (ft)  dp  (ft)  /ftp  exp  \j2n  (ft  +  nfi  +  pf2)  t  -  j(f>] 


Tl—  —  OO  p=  —  OO 
00  00 


m2  E  E  Jn  (ft)  dp_i  (ft)  /ftp  exp  [j2ir  (ft  +  nfi  +  pf2)  t  +  j<f>\ 


Tl—  —  OO  p=  —  00 

00  00 


m2  E  E  Jn  (ft)  dp+i  (ft)  /ftp  exp  [j2tt  (ft  +  nfi  +  pf2)  t  -  j<f>\ 


Tl—  —  OO  p=  —  OO 


This  simplifies  to  a  compact  expression  for  the  signal  after  the  filter  in  terms  of  its 
frequency  components, 


OO  OO 


^ det 


( t )  =  Kyj2PovtRe  E  E  ij,pexp  b'2^  (/c  +  n/i  +  p/2)  t] 


(2.6) 


n= — 00  p=— 00 
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where  I  define 


fn,r=K,P{Jn(MJAM  (2.7) 

+  \Jn- 1  {Pi)  +  Jn+l  (/dl)  6  Jp  (/?2 ) 

+  -m2  Jn  (/3i)  [Jp-i  (/?2)  +  </p+i  (/3i)  e 

The  electric  field  is  incident  upon  a  photodetector  at  the  termination  of  a  fiber-optic 
link.  The  photodetector  is  assumed  to  be  an  ideal  square-law  detector  operating  in 
its  linear  regime  with  responsivity  3?  .  The  photocurrent  is 


OO  OO  OO  OO 


iz{t)  =  ?RPopt  l  X  X  X  inJa%  exP  b'2yr  ([n-g\fi  +  [p~  k]  f2 )  t] 

n=— oo  p=—og  g=— oo  k=—oo 

(2.8) 

This  can  be  split  up  into  the  dc  term,  harmonics  of  fi,  harmonics  of  /2  and  mixtures 
between  f\  and  /2. 


OO  OO 


p(t)=sp,P  5]  5]  fcf 

L  g=—oo  k=— oo 
oo,n/g  oo  oo 

+  X  X  X  exP  b’27r  b  ^  0]  /i*] 

n=— oo  g=—oo  k=— oo 
oo,p^k  oo  oo 

+  X  X  X  exP  b'27r  [p  -  *] 

p=— oo  g=—oo  k=— oo 
oo,n^goo,p^k  oo  oo 

+  X  X  X  X  eXP  b'27r  (fa  -  s]  /l  +  fa  -  /2)  *] 

n= — oo  p=—oo  g=— oo  k=—oo 

The  indices  of  each  infinite  sum  can  be  renumbered  to  obtain 


OO  OO 


i‘(t)=!R/V  V  V  |j^|2 

^  g=—oo  k=— oo 
00,717^0  oo  oo 

+  E  E  E  jn+g,kfg*k  eXP  [j2^/^] 

n=—oo  g=— oo  fc=— oo 
oOjPt^O  oo  oo 

+  E  E  E  jg,p+kfg*k  eXP  b'27Tp/2t] 

p=— oo  g=— oo  k=—oo 

oo,n^0  oo,pt^O  oo  oo  ^ 

”  X  X  X  X  in+g,P+kjg*k  eXP  h'27r  fa/l  +  P./2)  t]  > 

n=—oop=—oog=—ook=—oo  ) 

The  double  infinite  sums  over  n  and  p  are  rewritten  as  singly  infinite  sums,  and  the 
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sums  over  negative  integers  have  their  signs  flipped  giving 


X  E  \fj2 

\9=— oo  k=—oo 
oo  oo  oo 

+  Un+g.kfg*k  exp  [j^nnfg]  +  jln+g,kjg*k  exp  [-j2'Knflt]) 

n=  1  g=— oo  k=— oo 
oo  oo  oo 

+  E  E  E  te,P+feig,fcexP  b'2^P/2^]  +  jg,-p+kfg*k  exp  [~j2npf2t]) 

p=  1  g=— oo  k=—oo 
oo  oo  oo  oo 

+  EE  E  E 

n=l  p=l  g=—oo  k=— oo 

{j n+g,p+kJ g*k  eXP  [/27r  («/l  +  Vf 2)  *]  +  j-n+g,-p+kjg*k  exp  [— j2tT  (n/i  +  pf2 )  t] 
+jn+g,-p+kfg*k  exp  b27r  (w/i  -  p/2)  *]  +  j-n+g,p+kfg*k  exp  [-/2tt  (n/i  -  p/2)  t])  } 

A  number  added  to  its  complex  conjugate  is  twice  the  real  part.  With  this  simpli¬ 
fication,  this  arranges  to  a  final  expression  for  the  photodetector  output  given  an 
arbitrary  filter: 


m  =SP^JJe  l  V  5]  jyt| 2  (2.9) 

V  g=—oc  k=— 00 
00  00  00 

+  25^  .?/+3:fcCfcexPb27rw/ii] 

7i=l  g=—oo  k— — oo 
00  00  00 

+  2EE  e  fg,p+kjg*k  exp  b2*P/2f] 

p=l  g=— 00  k=—oc 
00  00  00  00 

+  2EEE  e  jn+g,p+k3g*k  exp  b27r  (w/i  +  ph)  t] 

n= 1  p=l  g=—oo  k=—oo 

00  00  00  00  ^ 

+2  ^  jn+g,—p+kjg*k  exP  [j27r  (nfl  ~  Ph)  A  \ 

n=  1  p=l  g=—oc  k=— oo  J 

The  double-sum  over  indices  g  and  k  gives  the  contribution  of  each  pair  of  optical 
sidebands  that  beat  together  to  produce  the  rf  photocurrent.  In  this  form,  the  current 
is  separated  into  different  frequency  components  which  are  indicated  by  the  summa¬ 
tion  indices  n  and  p  .  The  first  term,  where  n  and  p  are  both  identically  zero,  gives  the 
dc.  The  second  term,  a  summation  over  the  index  n,  gives  the  fundamental  tone  at 
frequency  f\  and  its  harmonics.  The  third  term,  a  summation  over  the  index  p,  gives 
the  fundamental  tone  at  frequency  f2  and  its  harmonics.  The  fourth  term  is  the  sum 
frequencies  produced  by  the  mixing,  and  the  fifth  term  is  the  difference  frequencies 
produced  by  the  mixing. 
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2.4.2  Small  signal  approximation 

For  small  modulation  depth,  (3  <C  1,  and  no  residual  IM,  m  =  0,  the  Bessel  func¬ 
tions  can  be  approximated  by  Jq(/3 )  ~  1  and  Jn(/3)  ~  (/3/2)^  /  |n|!,  for  positive  n, 
noting  that  J-n(/3 )  =  (—l)nJn(/3).  Keeping  terms  of  lowest  polynomial  order,  the 
current  simplifies  to  the  following  equation  (2.10).  This  equation  gives  the  small 
signal  approximation  for  any  frequency: 


iz{t)  =MPoptRe\\hz0 


(2.10) 


00  n  Rn  ( _ 1  \9 

2  E  E  &  (n  -  g)\g\h”-g> ob-gfiVW 


n=l  g= i 
oo  P  op 


pi  (-i r 

2 ip  (p  -  fc) 

yy  yy _ p_  pnpp 


2  E  E  wb-kwK*-kK*-k  exp  [j27rp/2t] 

p=  1  k= 0  ^ 


(-1) 


^EEEE  2«+p  (n  —  )i  I  {  kyk]hn-g,P-kh-g,-k  exp  [j2vr  (nA  +  pf2)  t] 

n= 1  p=l  c/=0  fc=0  V  W 

(_-|\P+9+fc 


2=1  p=l  g=U  fc=i 

oo  oo  n  p  RnRV 

^  ^  2n+p  (n  —  g)\g\  (p  —  k)\k\ 

n=  1  P=1  5=0  fc=0  v  y  ) 


hzn_g_p+khzPgk  exp  \j 2tt  (n/i  -  pf2)  t] 


There  are  four  current  components  of  interest.  The  amplitude  of  the  dc,  as  should 
be  expected,  is  proportional  to  the  optical  power  in  the  optical  carrier  after  the 
filter.  The  current  at  the  fundamental  frequency  f±  is  linearly  proportional  to  the 
modulation  depth.  It  depends  on  the  negative  and  positive  first-order  sidebands 
beating  with  the  optical  carrier.  The  current  at  the  second-harmonic  frequency  2  f\ 
has  a  quadratic  relationship  to  modulation  depth.  It  depends  on  the  second-order 
sidebands  beating  with  the  optical  carrier,  as  well  as  the  first-order  sidebands  beating 
with  each  other.  The  current  produced  at  the  difference  frequency  2/i  —  f2  is  a 


third-order  intermodulation  product.  These  currents  are 

4  =»PoP,xs  (2.ii) 

t)1  {Xp  exp  [j2jrfit]}  (2.12) 

i2fl  =^Popt\plRe  {Xz  exp  [j4vr fR}}  (2.13) 

i2fl-h  =^Popt\plP2Re  {Xz  exp  [j2vr  (2 h  ~  A)  *]}  (2.14) 


where  for  convenience,  I  dehne  the  following  complex  constants,  which  I  will  call  link 
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distortion  constants: 


yz  _ uz  uz* 

~'Lo,oll'0,o 

(2.15) 

vz  _ ~L.Z  1  Z*  LZ  1  Z* 

A!  “/7%Er=-oot;+s,_p+fc4A°  0’°  °>°  -i>° 

(2.16) 

Yz  —hz  hz*  4-  hz  hz* 

(2.17) 

yz  _UZ  1  Z*  9;  Z  7  Z*  1  7  Z  7  z* 

^2  —U2,0U0,0  Z'il,0'i-1,0  “T  n 0,0'^— 2,0 

(2.18) 

yz  _  uz  UZ*  \  hz  hz*  _L_  9/7 2  hz* 

^-3  —  U2,-1U0,0  '  u2,0u0,l  "T" 

(2.19) 

_l_  hz  hz*  —  hz  hz*  —  2 hz  hz* 

'  't0,0'i-2,l  'i0,-l'i-2,0 

For  a  balanced  detector  system,  the  currents  subtract  from  each  other.  The  link 
constants  for  each  branch  can  be  subtracted  from  each  other  such  that  X0  =  X(J  —  X(2  , 
X\  =  X\  —  Xl,  etc.  Each  rf  photocurrent  outputs  an  rms  power,  which  is  proportional 
to  the  square  of  the  dc  current,  into  the  load  impedance,  Zout.  The  powers  for  the 
signal,  second  harmonic,  and  intermodulation  distortion  are  as  follows: 


Ph  =t|Z„,|»2P0^|V1|2 

(2.20) 

P2h  |Z„,|  |X2|2 

(2.21) 

Pih-h  \Zout\ ifOM  IWI2 

(2.22) 

In  this  section,  I  have  derived  closed  form  expressions  for  the  photocurrents  at 
different  frequencies  at  the  output  of  a  filtered  FM  link.  A  general  result  has  been 
given  in  (2.9)  which  includes  residual  intensity  modulation,  and  can  be  solved  to 
arbitrary  precision  by  taking  a  large  number  of  terms  in  the  infinite  sum.  A  small 
signal  approximation,  (2.10),  gives  the  output  current  at  any  frequency  component 
of  interest.  Expressions  for  the  photocurrent  at  the  fundamental,  second  harmonic 
and  third  order  intermodulation  distortion  have  been  derived,  which  will  be  useful  in 
expressing  figures  of  merit  for  distortion  and  dynamic  range. 


2.4.3  Gain 


For  an  PM  link,  the  modulator  is  driven  by  an  applied  voltage.  The  peak  input 
voltage,  Vin  produces  an  rms  input  power  Pin  when  delivered  to  a  load  impedance  Zin 
such  that  f3\  =  (-7T V)n/14-)2  =  7t22 Pin  \Zin\  /V%.  The  output  signal  is  given  by  (2.20). 
The  gain  is  therefore 


G 


PM  ~ 


(2.23) 


For  an  FM  link,  the  modulation  efficiency  of  a  current  modulated  FM  laser  is  r/,  in 
units  of  Hz/A,  typically  of  the  order  of  a  few  hundred  MHz  per  m  A.  The  peak  input 
current,  im  gives  a  modulation  of  <5/  =  (r]iin)  2  =  2 ifPin/  \Zin\.  The  gain  is  therefore 


I  7 

_  |  ^  out 

FM  ~  I  ~  I 
|  ^in  | 


(2.24) 
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2.4.4  RF  noise  figure 

In  this  section,  I  derive  the  signal  to  noise  ratio  (SNR)  for  the  small  signal  approxima¬ 
tion  of  an  arbitrary  link  and  the  noise  figure.  A  passive  link  with  no  amplification  will 
be  considered,  so  the  primary  noises  seen  at  the  detector  are  shot,  thermal,  phase  and 
RIN.  The  shot  noise  spectral  density  is  proportional  to  q,  the  elementary  charge,  and 
to  the  total  dc  from  the  photodetectors,  idc  =  3? Popt  l-^-o  !■  The  thermal  noise  spec¬ 
tral  density  is  equal  to  the  product  Boltzmann’s  constant,  kp,  and  the  temperature, 

Tk. 


Ssn  =2qidc\Zout\  (2.25) 

Stn  —kpJ^K  (2.26) 


Assuming  a  Lorentzian  model  for  the  laser’s  spectral  line,  the  phase  noise  on  the 
optical  carrier  is  white  noise  with  spectral  density  proportional  to  the  laser’s  3-dB 
linewidth,  An  [62],  The  phase  fluctuations  are  converted  to  intensity  fluctuations 
by  the  filter  in  the  same  manner  as  it  converts  the  modulation.  The  average  phase 
fluctuations  in  a  small  bandwidth  near  some  frequency,  /,  are 


(p  (tf) 


Az/A / 

7T  f2 


Near  the  first  modulation  frequency,  fi,  the  power  spectral  density  of  the  phase  noise 

is 

Spn  «  \Zoui\  M2PL^  |Ad|2  (2.27) 

TT/l 

The  modulation  is  assumed  to  be  below  the  relaxation  frequency  of  the  laser,  so  the 
RIN  is  modeled  as  white  noise.  The  power  spectral  density  of  the  noise  at  the  output, 
near  the  modulation  frequency  is 

\\ZM\WPppp-\Ytf  (2.28) 

where  B  is  the  bandwidth  in  Hz.  The  total  noise  power  is 


(2.29) 


The  noise  figure  is  given  by  relation  to  the  gain,  given  in  (2.23)  and  (2.24)  ,  and  for 
a  thermal  noise  limited  input  as 


NF 


1  (Ssn  +  Spn  +  Sin ) 
+  G  +  GkBTK 


(2.30) 


For  large  positive  link  gain,  approximations  to  the  noise  figure  expressions  for  shot 
and  phase  noise  limiting  are  given  in  Table  2.1  on  page  25. 

It  is  important  to  warn  that  the  effect  of  noise  caused  by  optical  amplification  has 
not  been  explicitly  included  in  this  analysis.  It  is  likely  that  optical  amplification, 
such  as  an  EDFA,  will  be  used  in  high  performance  PM/FM-DD  link  architectures 
since  the  NF  and  spurious  free  dynamic  range  scale  with  the  optical  power.  The  noise 
degradation  from  the  laser’s  intrinsic  noise  by  the  amplifiers  must  be  included  in  the 
quantity  provided  for  the  laser  phase  and  amplitude  noise. 
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Table  2.1:  Approximations  to  the  noise  figure  expressions  for  arbitrarily  filtered  links. 
These  assume  large  positive  gain  with  either  shot  or  phase  noise  limited  noise  figures. 
Shot  noise  limits  occur  for  moderate  optical  powers  and  phase  noise  limit  occurs  for 
much  larger  optical  powers.  These  approximations  are  not  valid  if  the  link  attenuates 
the  rf  power. 


PM 

FM 

Shot  noise  SNR 

MPoPt/3  fIXil2 

MPoptS2fl  |Xi|2 

4gB£  X§ 

4 qBfl  £  X§ 

Phase  noise  SNR 

"SIR 

2A  uB 

nsh 

2A  uB 

Gain 

1  7  I  I  ry  \  (  7 T^iPopt  1  y  \  \ 

l^inl  l^out]  l  y ^  1^-1 1  I 

\Zout\  { V^Popt  I  V'  1^ 

\Zin\  V  h  1  ) 

12 

1 

Shot  noise  NF 

1  ,  2qV*Z 

~T“  \Zin\fi2qT,\x£ 

\Zin\ytP0ptir2\X1\2kBTK 

+  3tPovtV2\Xi\2  kBT 

K 

Phase  noise  NF 

“r  \zinW3RkBTK 

1  |  \Zin\^V 

ri2irkBTK 

2.4.5  Distortion 

The  signal  distortion  caused  by  the  link  can  be  described  by  the  output  power  at 
frequencies  that  are  harmonics  and  mixing  terms  of  the  modulation  frequencies.  For 
now,  I  assume  there  is  no  residual  IM,  and  assume  two  modulation  tones  have  equal 
modulation  depth,  f3  —  /3\  —  fa  for  PM,  or  sf  =  Sh  =  Sf2  for  FM.  The  second-order 
output  intercept  point  (OIP2)  and  third-order  output  intercept  point  (OIP3)  for  PM 
and  FM  are  calculated  in  the  table  below. 


2.4.6  Spurious  free  dynamic  range 


The  spurious  free  dynamic  range  (SFDR)  is  defined  as  the  SNR  at  the  maximum 
usable  modulation  depth.  This  can  be  defined  when  either  the  second-order  or  third- 
order  distortion  products  breach  the  noise  floor.  The  second  harmonic  is  equal  to 


the  shot  noise  power  at  modulation  depth  =  prq 


gBJ2\x§ 

5R  Popt 


and  to  the  phase 


noise  power  at  modulation  depth  (3f  = 


2Ai /B 


the  shot  noise  power  at  modulation  depth  = 


|Ad|//i|X2|. 

/256qBE|Vg| 


SRPe.pt  1-X  3r 

2  _  (  128Ayg|Xi|2  ^  1//3 


The  IMD3  is  equal  to 

1/3 

and  to  the  phase 

noise  noise  power  at  modulation  depth  /3f  =  J  .  Using  the  previously 

calculated  expressions  for  the  SNR,  the  spurious  free  dynamic  ranges  are  compiled 
in  Table  2.2.  In  addition,  the  third-order  limited  SFDR  could  also  be  calculated 
given  the  noise  figure,  gain,  and  output  intercept  points  in  dB  units  by  the  expres¬ 
sion  SFDR  =  |  (OIP3  -  G  +  174  ^  -  10  log10  (B)  -  NF ) .  These  figures-of-merit 
are  often  defined  with  respect  to  1  Hz  bandwidth.  They  generally  depend  on  the 
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Table  2.2:  General  expressions  for  01  P 2 ,  0IP3,  and  spurious  free  dynamic  range 
for  an  abitrarily  filtered  link  with  either  phase  or  frequency  modulation  and  direct 
detection  given  in  terms  of  the  link  distortion  constants.  SFDR  is  limited  by  either 
shot  or  phase  noise,  and  second-order  or  third-order  distortion. 


particular  modulation  frequencies  chosen.  Maximizing  the  ratios  of  X]  |  /  | X2  and 
|Ad|  /  |X3|  will  improve  the  dynamic  range  of  the  link. 

2.5  Mach-Zehnder  interferometer 

The  simplest  filter  used  as  a  discriminator  is  an  asymmetrical  Mach-Zehnder  inter¬ 
ferometer  (a-MZI).  I  derive  figures  of  merit  for  this  link  in  order  to  verify  the  general 
theory  against  previously  published  results.  One  arm  of  the  interferometer  has  a  time 
shift  with  respect  to  the  second  arm,  denoted  by  r  .  We  assume  50%  coupling  ratios 
and  quadrature  bias,  obtained  by  choosing  the  carrier  frequency  and  time  delay. 

The  filter  transfer  functions  seen  by  the  two  output  branches  of  the  Mach  Zehnder 
interferometer  are 

K,P  =  \~J2exp  [_j27r  (n-ft  +  P &  (2-31) 

K,P  =  |  \  exp  [~j2n  (nfi  +  pf2 )  r]  (2.32) 

By  taking  the  absolute  value  squared  of  either  transfer  function,  one  can  see  that  its 
intensity  response  is  the  familiar  raised  sinusoid  and  it  is  quadrature  biased  with  half 
the  carrier  power  transmitted  to  each  branch. 

Using  the  transfer  functions,  I  evaluate  the  link  distortion  constants  for  both 
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branches: 


Xl  =- 
0  2 

(2.33) 

xi  =\j  (!  -  exp  [-j27t/it]) 

(2.34) 

xi  =\  (1  +  exp  [-j2nf1r]) 

(2.35) 

X]  =0 

xl  =  -  4  sin2  [vr/ir]  sin  [nf2r]  exp  [-ji r  (2/i  -  f2)  r) 

(2.36) 

(2.37) 

and 

xo  =~ 

(2.38) 

X'f  =  -  (1  -  exp  [-j2TTf1r\) 

(2.39) 

if  (1  +  exp  [-j27t/it]) 

(2.40) 

Af  =0 

x3  =4  Sin2  [vr/ir]  sin  [7 r/2r]  exp  [-jn  {2fx  -  f2)  r] 

(2.41) 

(2.42) 

As  expected  for  an  MZI  at  quadrature,  I  find  that  there  is  no  second-harmonic  so 
that  OIP2  is  infinite.  For  the  FM  link,  we  choose  a  short  time  delay  such  that 
approximation  f\T,  f2r  <C  1  is  valid.  The  absolute  value  of  the  other  coefficients  after 
the  balanced  detection  are 

Ah  =2  sin  (vr/ir)  ~  27t/i r 

|A3|  =8 sin2  (vr/ir)  |sin(7r/2r)|  «  8n3tff2T3 

(2.43) 

(2.44) 

A  summary  of  the  hgures  of  merit  are  given  in  the  table  below.  The  same  results  are 
found  by  [17],  which  supports  the  general  analysis.  The  important  result  from  [17] 
was  that  the  shot  noise  limited  spurious  free  dynamic  range  of  the  PM-MZI  link  is 
identical  to  that  of  a  Mach  Zehnder  modulated  IM-DD  link. 

2.6  Complementary  linear-field  demodulation 

In  this  section,  I  discuss  filter  transfer  functions  that  allow  for  highly  linear  dis¬ 
crimination.  I  find  that  the  ideal  system  has  two  filters  with  ramps  of  electric  held 
transmission  versus  frequency,  and  linear  phase. 

A  number  of  groups  have  proposed  or  built  optical  filters  that  have  a  transfer 
function  linear  in  optical  intensity  versus  frequency  and  small  group  delay.  Within 
one-half  period,  the  transfer  function  can  be  represented  by 

K,P  =  VA(fb  +  nfi  +  p/2)  exp  [— j'2tt  (fb  +  nf 1  +  pf2 )  r }  (2.45) 
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Table  2.3:  Figures  of  merit  for  an  PM-DD  link  with  an  a-MZI  and  balanced  detection. 


PM 

FM  (small  delay) 

Gain 

\Zin\  \Zout\ 4  sin  (7t/it) \j 

fcl 

Shot  noise  NF 

1  1  qV« 

1  |  1  Zin | Q 

\Zin\2idcn2kBTK  |sin(7r/ir)k 

_l"  3tPoptV22ir4T2kBTi{ 

Phase  noise  NF 

~ n  Avv- 

\Zin\ir3f{kBTK 

1  | 

ri2irkBTK 

OIP3 

A  7  1  ,-2  sin(7r/ir) 

1  out  \  f'dc  sin(7r/2r) 

4  -Zout 

Shot  noise  SFDR3 

f  sin4(7r/ir)2idc\2/3 
\  |sin(7r/2r)|qB  I 

^/fr2  (|f)  2,3 

Phase  noise  SFDR3 

(  ^  1 
^  Ai/S|sm(7r/ir)sm(7i'/2r)|  ) 

2/3 

1 

/  1  \  2/3 

\  AisBttt2  ) 

where  A  is  a  slope  in  units  of  inverse  frequency,  fb  is  a  bias  frequency  offset  from  the 
carrier,  and  r  is  a  time  delay  giving  the  filter  linear  phase.  The  intensity  response  is 


hn,PK,P  =  A  (fb  +  nf 1  +  pf-2 )  (2.46) 

which  is  linear  in  slope  A.  Using  the  transfer  function,  I  evaluate  the  link  constants: 

=y/Afbe-j2nfbT  (2.47) 


I  (249) 

-2f+lf~¥)  e~si’hT 

x> =Ah  <2’5o) 

-  J^j  e~P^h-h)  (2.51) 


fl 


—  e 


-j2nf!T 


fb 


(2.48) 


Ad  =Afb 


Generally,  X2  and  X3  are  non-zero  for  this  discriminator,  even  if  the  square  roots 
are  expanded.  This  means  that  a  discriminator  that  is  linear  in  optical  intensity  will 
still  produce  second-order  and  third-order  distortion.  Mixing  in  the  photodetector 


produces  cross  terms  that  are  not  eliminated.  An  FM  discriminator  that  is  linear  is 
optical  intensity  will  not  produce  a  distortion-less  link. 

The  ideal  discriminator  for  the  link  is  a  pair  of  optical  filters  that  are  linear  in 
electric  field.  Within  one  period,  the  field  transmission  ramps  linearly  with  frequency, 
and  the  filter  has  linear  phase.  The  transfer  functions  near  the  carrier  are 

K,p  =  A  ( fb  +  nfi  +  pf-2 )  exp  [—j2n  (fb  +  nfi  +  pf2)  r\  (2.52) 

hl,p  =  (fb  ~  nfi  -  pf2 )  exp  [~j2tt  (fb  +  nfi  +  pf2 )  r]  (2.53) 

where  A  is  a  slope  in  units  of  inverse  frequency  and  r  is  a  time  delay.  The  1/V2 
prefactor  is  an  optical  splitter  before  two  physical  filters.  I  define  the  constant  T  to 
describe  the  dc  bias  of  the  filter,  which  is  the  fraction  of  optical  power  transmitted 
by  the  filter  at  the  optical  carrier  frequency.  The  link  distortion  constants  are 


*0  =A2fl/2  =  T/2 

(2.54) 

X\  =AfiT1/2e~j2nflT 

(2.55) 

(2.56) 

=A2f2e~j4nflT 

(2.57) 

xl  =o 

(2.58) 

X2  —T/2 

(2.59) 

X2  =  -  AfiTl/2e~j2nflT 

(2.60) 

y '2  —Ye-o^hT 

(2.61) 

X2  =A2  fle~jAirflT 

(2.62) 

xl  =0 

(2.63) 

All  higher  order  link-constants  are  zero.  The  non-zero  values  of  X/  are  due  to  the 
squaring  of  the  signal  at  the  detector.  The  distortion  is  caused  by  the  first-order  side¬ 
bands  beating  with  each  other.  However,  because  the  second  harmonics  are  in  phase, 
they  cancel  at  the  balanced  detector,  giving  perfect  distortionless  performance.  The 
current  component  at  the  fundamental  frequency  will  be  180°  out  of  phase  between 
the  two  photodetectors,  but  the  second-harmonic  will  be  in  phase.  After  the  balanced 
detector,  the  only  term  that  does  not  cancel  is 

\Xi\  =2AfiTl/2.  (2.64) 

It  is  important  to  note  that  the  intensity  modulation  term  also  cancels  because  of 
balanced  detection.  Residual  intensity  modulation  of  the  laser  and  relative  intensity 
noise  present  before  the  demodulation  will  not  be  present  at  the  output  of  this  system. 

In  the  small  modulation  depth  approximation,  this  ideal  link  has  no  other  higher- 
order  distortion.  Using  a  symbolic  algebra  solver,  I  verified  that  the  current  is  zero  for 
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all  intermodulation  and  harmonic  frequencies  up  to  sixth  order.  At  a  given  harmonic, 
sum  or  difference  frequency,  if  all  the  sidebands  in  the  sum  in  (2.9)  corresponding  to 
that  frequency  fall  within  a  region  of  the  filter  that  closely  approximates  the  desired 
linear  ramp  function,  the  output  current  is  zero. 

Additional  sources  of  nonlinearity  are  the  frequency  modulated  laser  source,  opti¬ 
cal  fibers  and  photodetector.  For  sufficient  modulation  depth,  the  dominant  sidebands 
will  fall  outside  the  bandwidth  of  the  filter  and  this  saturation  will  cause  nonlineari¬ 
ties. 

2.6.1  Noise  and  gain 

In  this  section,  I  will  consider  the  effect  of  the  bias,  T,  on  the  noise  figure  of  the  link. 
Low  biasing  the  filter,  to  decrease  the  dc  current  at  the  detector,  had  been  suggested 
by  [54]  and  others  to  improve  the  noise  figure  (NF)  of  a  PM  or  FM  link.  However, 
there  is  a  tradeoff  between  decreasing  the  dc,  which  decreases  shot  noise,  and  reducing 
the  signal  gain,  so  an  optimal  bias  point  must  be  found.  The  filter  cannot  be  biased 
exactly  at  the  null  or  the  link  would  have  zero  output  current,  since  I  find  in  (2.64) 
that  the  output  is  proportional  to  the  square  root  of  the  bias.  This  is  consistent  with 
experience  with  carrier  suppression  on  IM-DD  links. 

The  noise  figure  of  the  link  is  comprised  of  a  term  for  an  attenuated  link,  the  shot 
noise  component,  and  the  phase  noise  component.  Intensity  noise  does  not  appear 
because  it  is  canceled  with  the  balanced  detection.  The  noise  figures  for  PM  and 
FM  are  given  by  as  follows.  They  are  written  in  terms  of  the  dc  photocurrent  at  the 
detectors,  instead  of  the  total  optical  power  before  the  filters,  since  current  handling 
of  the  diodes  is  usually  a  limiting  factor. 

TV 2 

N FPM  =1  +  —  — 1/2-2  a2  n 

\Zin\  Fouf  47T 

<V»r  a v_yl 

\zin\ tdcP-2kBTKA^f;  \Zin\ flkBTK 
\Zin\  T 

nffm  =1  +  ■„  '|Z'9.2  A2 
\Zout\  4 ncA2 

\zin\  gT  \Zin\  Au 

idcri22kBTKA 2  rf2nkRTK 

A  useful  question  is  whether  it  makes  sense  to  low  bias  the  filter  in  an  attempt 
to  improve  the  noise  figure.  The  answer  depends  on  whether  the  designer  is  limited 
by  optical  power  available  or  by  the  maximum  photocurrent  the  photodetectors  can 
handle.  For  a  fixed  current,  for  which  the  optical  power  is  increased  to  maintain,  the 
derivative  of  the  NF  with  respect  to  the  bias  is 

dNFPM  _  V2  qV2 

d T  \Zin\  \Zout\ An2i2dcA2 f2  +  \Zin\  idcn22kBTKA2f2 

DNFfm  _  | Zin |  \Zin\q 

d  T  \Zout\Ar]2i2dcA2  idcq22kBTKA2 


(2.67) 

(2.68) 


(2.65) 


(2.66) 
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Table  2.4:  Gain  and  noise  figure  expressions  for  the  complementary  linear-field  de¬ 
mo  dulated  PM-DD  link. _ 


PM 

FM 

Gain 

\Zin\  \Zout\  4T_1  ( ^Af,y 

T-'  (^Af 

Shot  noise  NF 

1  1  qV* 

1  1  \Zin\q 

^  \Zin\^Popt-K22kBTKA2f2 

1  'r  SR PoptTl22kBTKA 2 

Phase  noise  NF 

i  |  AiAtf 

~l~  1  Zin\n3f?kBTK 

1  |  \Zin\^V 

ri2irkBTK 

Laser  linewidth  limit 

(NFpn  -  1)  1  Z^fPBTK 

_ Lze _ 

WFm  ~  1)  ,£tSt 

These  are  always  positive  quantities  so  the  noise  figure  monotonically  decreases  with 
decreasing  bias  as  long  as  the  current  is  maintained.  However,  the  phase  noise  will 
begin  to  dominate  over  the  shot  noise  when  Tq  <  2 idcA2^-,  and  any  NF  improvement 
will  be  negligible.  For  example,  with  Av  =  1  MHz  and  A  =  1/50  GHz  ,  choosing  a 
low  bias  point  only  makes  sense  if  the  maximum  dc  current  is  less  than  160  /i A. 

If  the  available  optical  power  is  fixed,  idc  =  3? PoptT ,  then  the  derivative  of  the  NF 
with  respect  to  the  bias  is  always  negative: 

ONFpm  _ _ Vi _ 

dT  \Zin\  |  Zout  4n2WP2ptT2A2f? 

3NFfm  _  |  Zin  | 

8T  ~  ~  \  Zout  \  4q2^R.2  P2ptT2  A2 

Reducing  the  bias  only  serves  to  reduce  the  gain  of  the  link  and  the  NF  gets  worse 
with  the  lower  bias. 

For  high  optical  powers,  the  NF  is  phase  noise  limited.  This  is  independent  of 
the  filter  bias  and  the  slope  of  the  filter.  Because  random  frequency  fluctuations 
are  added  to  the  optical  carrier  at  the  same  time  as  it  is  modulated,  the  maximally 
achievable  SNR  is  set  at  the  laser,  and  cannot  be  improved  by  the  rest  of  the  system. 
This  formula  sets  a  fundamental  relationship  between  the  maximally  achievable  noise 
figure,  the  laser  linewidth  and  the  modulation  efficiency.  For  a  given  noise  figure  and 
modulation  efficiency,  the  maximum  laser  linewidth  is  given  in  the  table  above.  This 
fundamental  relationship  between  modulation  efficiency,  linewidth  and  noise  figure  is 
plotted  in  Fig.  2.5  for  a  range  of  typical  values. 

2.6.2  Transfer  function  curvature 

The  realized  complementary,  linear-field  filters  will  be  non-ideal  in  their  phase  and 
amplitude  transfer  functions.  It  is  useful  to  define  a  mask  or  envelope  for  the  filter 
transfer  function’s  phase  and  amplitude  which  will  guarantee  a  desired  nonlinearity 
performance.  I  find  that  a  mask  with  the  form  of  a  quadratic  equation  of  frequency 
offset  from  the  carrier  will  give  a  tractable  bound  for  OIP2,  and  a  cubic  mask  will 
give  a  tractable  bound  for  OIP3 .  Local  derivatives  do  not  affect  distortion  as  long  as 
the  function  falls  within  given  envelope. 


(2.69) 

(2.70) 
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Figure  2.5:  Phase  noise  limited  noise  figure  versus  linewidth  and  modulation  efficiency, 
assuming  a  50  ohm  impedance. 

Phase  noise  limited  noise  figure  (dB) 


The  realized  transfer  function  for  one  branch  of  the  discriminator  is  written  in  the 
form 

Kf)  =  ^  (Vf  +  Af  +  Aa  (/))  exp  [~j2nfr  -  jAp  (/)]  (2.71) 

where  Aa  (/)  and  Ap  (/)  are  the  deviations  from  the  ideal  phase  and  amplitude,  and 
/  is  the  offset  from  the  carrier. 

Figure  2.6  on  page  33  illustrates  the  masks  for  the  amplitude  and  phase  for  bound¬ 
ing  the  second-order  figures  of  merit.  The  deviations  from  ideal  for  the  amplitude  and 
phase  must  fall  within  bounds  which  relax  further  away  from  the  carrier  frequency: 


and 


Aa  (/)  =  £2  (/)  A2/2,  (2.72) 

M/)l  <  e2,max  (2.73) 

A  P(/)=<M2/2  (2.74) 

|02  (/)  |  <02, max  (2.75) 


where  £2, max  and  02, max  are  small  positive  constants.  For  a  two-tone  test  derivation,  I 
make  the  approximations  that  the  modulation  tones  are  closely  spaced,  f\  ~  ,/2  =  /, 
the  phase  deviation  is  small  so  that  exp  \—j(f>]  ~  1  —  j(j)  and  the  frequency  is  low  with 
respect  to  the  bias  so  that  Af  <  y/T.  For  the  OIP2  derivation,  1  use  the  second 
harmonic  as  the  distortion  term  of  interest.  1  also  assume  complementary  filters  and 
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Figure  2.6:  Illustration  of  the  quadratic  envelope  on  the  transfer  function  that  bounds 
the  second-order  figures  of  merit  for  the  complementary  linear-field  discriminator. 

Quadratic  envelope  of  field  amplitude  Quadratic  envelope  of  phase 


balanced  detection.  The  two  tone  transfer  function  for  one  branch  is 

hln,P  =^=  (yT  +  A  (nfi  +  pf2 )  +  e2  (nfi  +  pf2 )  A2  (nf\  +  pf2)2 )  (2.76) 

•  exp  [-j2vr  (nfi  +  pf2)  r  -  jf2  (nfi  +  pf2 )  A2  (nfi  +  pf2 )2] 

After  algebraic  simplihcations,  assuming  the  worst  case  addition  of  errors,  the 
second-order  link  distortion  constant  is  bounded  by 


where 


|X2|  <  A2 f 2  (Cie2  ,max  +  jc2<h  ,max) 

Ci  =12Vf 


2  r  2 


C2  =12  T  -  4  A2f 


and  its  magnitude  is  therefore 


\X2\  <  eimaxT  + 


T  -  -  1-  f! 


The  second-order  output  intercept  point  is  lower  bounded  as 

8  1 

OIP2  >  -Rloadtic - 7 

~9  elmaxT  +  flmax(T-\A*ry 

for  the  worst  case  frequency, 

OIP 2  >  jj Rioadil  1 


cf2  T  A-  (h2  T2 

fc2,mnii  '  <r2,max± 


(2.77) 

(2.78) 

(2.79) 


(2.80) 


(2.81) 


(2.82) 
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Figure  2.7:  Illustration  of  the  cubic  envelope  on  the  transfer  function  that  bounds 
the  third-order  figures  of  merit  for  the  complementary  linear-field  discriminator. 

Cubic  envelope  of  field  amplitude  Cubic  envelope  of  phase 


and  for  the  important  half  field-bias  case  where  y/T  =  '/a, 

32  1 

OIP2  >  - Rloa, *il-2 - -T-y -  (2.83) 

^  ^2, max  ^2, max 

Figure  2.7  on  page  34  illustrates  the  masks  for  the  amplitude  and  phase  for  bound¬ 
ing  the  third-order  figures  of  merit.  The  deviations  from  ideal  for  the  amplitude  and 
phase  must  fall  within  bounds  which  relax  further  away  from  the  carrier  frequency: 


and 


Aa  (/)  =  £3  (/)  A3f3,  (2.84) 

M/)  \<e3,max  (2.85) 

Ap  (/)  =03  (/)  A3f3  (2.86) 

1 03  (0)  |  <03 ,«  (2.87) 


where  £3, max  and  4>3:max  are  small  positive  constants.  For  a  two-tone  test  derivation, 
I  make  the  same  approximations  as  before.  The  two  tone  transfer  function  for  one 
branch  is 

=~^  (yT  +  A  (n/i  +  pf2 )  +  £3  (nfi  +  p/2)  7l3  (n/i  +  p/2)3) 

•  exp  [— j'2tt  (nfi  +  pf2 )  r  -  j03  (n/i  +  pf2)  A3  (n/i  +  p/2)3]  (2.88) 

After  algebraic  simplifications,  assuming  the  worst  case  addition  of  errors,  the 
third-order  link  distortion  constant,  given  by  is  bounded  by 
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where 


X3  <  A3f  (Cl£ 3  ,max  +  jC 2^3  ,max) 


(2.89) 


Ci  =24  Vf 
C2  =24  T  +  36  A2  f2 


and  its  magnitude  is  therefore 


|A'3|  <  A3  f  24,^  +  <g~T  (l  +  ?A2/2/t) 

The  third-order  output  intercept  point  is  lower  bounded  as 

12  i 


OIP3  >  rRioadf 


for  the  worst  case  frequency, 


+ AAA  (i + §.42/77) : 


4  x 

OIP3  >  -Rload  m 

-  3  T 


2 

dc 


yj £l,max  +  ^3 ,maxT  (1  +  3/32T)" 


and  for  the  important  half  field-bias  case  where  \/T  =  x/2, 

16  1 

01  P3  >  —  R-loaSdc  I 

3  /p  2  I  ^2  121 

A/  a3,max  ~  ^3, max  256 


(2.90) 

(2.91) 


(2.92) 


(2.93) 


(2.94) 


(2.95) 


1  performed  Monte  Carlo  simulations  to  verify  these  error  bounds.  1  created 
a  complementary  linear-field  transfer  function  and  added  random  deviations  that 
fall  within  the  mask.  The  transfer  function  was  used  to  analytically  calculate  the 
distortion  figures  of  merit.  This  was  repeated  1000  times  for  each  parameter,  and 
the  worst  case  was  saved.  The  worst-case  simulated  distortions  fell  within  0.5  to  2 
dB  above  the  lower  bound,  making  this  a  suitable  mask.  The  best  cases  sometimes 
outperformed  the  bound  by  10s  of  decibels,  but  this  was  highly  dependent  on  the 
modulation  frequency. 


2.6.3  Residual  intensity  modulation 

Residual  intensity  modulation  sets  a  lower  limit  on  the  distortion  for  a  link  using 
complementary  linear  held  discriminators.  The  effect  of  residual  IM  can  be  obtained 
from  (2.9).  It  is  difficult  to  write  a  general  expression,  but  it  is  possible  to  expand 
some  individual  terms.  In  lowest  polynomial  order  of  the  modulation  depth,  the 
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Tabic  2.5:  Expressions  for  the  worst  case  OIP2 ,  OIP3,  and  spurious  free  dynamic 
range  for  complimentary  linear-field  demodulation  limited  by  filter  curvature. 
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Figure  2.8:  Monte  Carlo  simulation  results  to  test  the  suitability  of  the  derived  bounds 
on  the  0IP2  and  0IP3.  Each  point  is  the  worst  case  of  1000  trials  with  random  errors, 
and  is  compared  to  the  analytical  bounds.  We  assume  closely  spaced  tones  around  2 
GHz,  1/10  GHz  slope,  5  mA  of  current  per  detector  (fdc  =  10mA),  50  ohm  impedance, 
and  0.5  amplitude  bias,  T  =  0.25 .  The  analytical  expression  bounds  the  simulation 


within  less  than  2  dB. 
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currents  of  interest  are 


>/,  ^Pv,Re  {  [ft  (hlX%  -  *5,0*",  „) 
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(2.96) 

(2.97) 


(2.98) 

(2.99) 
(2.100) 


It  is  useful  to  normalize  the  IM  to  the  FM.  One  method  of  normalization  is  to 
look  at  the  optical  power  the  IM  and  angle  modulation  contribute  to  the  first-order 
optical  sidebands  in  the  small  signal  approximation.  The  optical  power  in  the  first 
order  sidebands  from  the  IM  is  Popt'm\/ 16.  The  optical  power  in  the  first  order 
sidebands  clue  to  the  PM  or  FM  is  Poptf32/4:.  I  define  a  relative  residual  IM,  T,  as 


r  = 


m 

2/3 


(2.101) 


The  corrected  expressions  for  the  link  distortion  constants,  including  the  effect  of 
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residual  intensity  modulation,  are 


X*  =hz00hz0*0  (2.102) 

Ai  =h10h00  —  h00h_10  (2.103) 
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For  the  complementary,  linear-field  demodulation,  the  magnitude  of  the  distortion 
constants  are 


|Xi|  =2Af1T1/2 

(2.106) 

\X2\  =8AfiT1//2T  cos  (0/2)| 

(2.107) 

A3  =4XT2T1/2  \2fi  +  /2exp  [j2(f)\\ 

(2.108) 

Since  the  intensity  modulation  is  residual,  the  frequency  modulation  will  be  much 
greater  than  the  intensity  modulation.  With  balanced  detection,  both  the  dominant 
second-harmonic  terms  and  dominant  IMD3  terms  are  quadratic  with  the  intensity 
modulation  are  linear  in  the  IM.  The  values  for  residual  intensity  modulation  limited 
OIP2  and  OIP3  are  in  the  below  table.  A  set  of  example  curves  are  shown  in  2.9.  It  is 
interesting  to  note  that  the  values  for  the  spurious  free  dynamic  range  are  independent 
of  the  bias. 

2.6.4  Dispersion 

The  dispersion  of  the  optical  fiber  also  increases  the  distortion  of  a  PM-DD  or  FM- 
DD  link.  The  dispersion  is  modeled  by  multiplying  the  filter  transfer  function  by  the 
term  exp[— jnDz  (nfi  +  p/2)2]  ,  where  D  is  the  fiber  dispersion  parameter  and  z  is 
the  fiber  length.  The  figure  below,  2.10,  shows  example  curves  of  the  upper  limit 
the  dispersion  sets  on  OIP3 .  It  degrades  by  20  dB  per  decade  of  fiber  length.  This 
can  be  corrected  by  using  a  length  of  dispersion  compensated  fiber,  or  by  designing 
a  discriminator  filter’s  transfer  function  to  include  the  inverse  of  the  dispersion.  The 
mechanism  for  the  dispersion’s  impact  on  the  link  distortion  is  conversion  of  phase  or 
frequency  modulation  to  intensity  modulation. 
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Tabic  2.6:  Expressions  for  01 P 2  and  0IP3  for  complimentary  linear-field  demod¬ 
ulation  limited  by  residual  intensity  modulation,  with  an  arbitrary  phase  difference 
between  the  angle  modulation  and  the  intensity  modulation.  The  frequency  depen¬ 
dent  terms  are  only  a  small  correction  for  closely  spaced  tones. 
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Figure  2.9:  OIP3  and  SFDRg  for  an  ideal  discriminator  for  different  values  of  residual 
intensity  modulation,  assuming  closely  spaced  tones  around  2GHz,  5  mA  of  current 
per  detector  (idc  =  10mA),  50  ohm  impedance,  and  0.5  amplitude  bias,  T  =  0.25. 
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Figure  2.10:  OIP3  for  complementary  linear-field  discriminators  for  different  slope 
values  and  fiber  dispersion,  assuming  standard  SMF,  with  D  =  — 20Ps2/km,  closely 
spaced  tones  around  2GHz,  5  mA  of  current  per  detector  ( idc  =  10mA),  50  ohm 
impedance,  and  0.5  amplitude  bias,  T  =  0.25. 
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2.7  Summary 

In  this  chapter,  I  have  proven  theoretically  that  complementary  linear-field  discrim¬ 
inators,  if  implementable  with  real  optical  filters,  can  potentially  lead  to  microwave 
photonic  links  with  very  high  dynamic  range.  Table  2.4  summarizes  the  noise  figure 
metric  in  the  shot  noise  and  phase  noise  limited  regimes,  table  2.5  gives  limits  on  the 
spurious  free  dynamic  range  by  filter  curvature,  and  table  2.6  gives  the  SFDR  lim¬ 
ited  by  residual  intensity  modulation.  Assuming  the  link  is  limited  by  photodetector 
current  rather  than  optical  power,  I  find  that  the  gain  and  noise  figure  both  benefit 
from  low  biasing  the  discriminators.  In  the  next  chapter,  the  arbitrary  filter  model 
derived  here  will  be  used  to  evaluate  physical  implementions  of  the  discriminators,  to 
predict  the  limits  of  their  performance. 
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Chapter  3 

Simulated  filter  performance 


Complementary  linear-field  demodulation  can  achieve  high  dynamic  range  if  good 
approximations  to  the  desired  filter  transfer  functions  can  be  physically  realized.  In 
recent  years,  there  has  been  much  work  in  devising  microwave  photonic  filters  [8,  63]. 
As  reviewed  by  [64] ,  a  systematic  way  that  microwave  photonic  filters  can  be  designed 
is  by  using  techniques  borrowed  from  the  field  of  digital  filters.  One  specifies  the 
coefficients  of  the  z-transform  representation  of  the  filter,  and  then  uses  a  synthesis 
algorithm  to  map  to  optical  components  such  as  couplers,  resonators,  and  delay  lines. 
The  problem  of  discriminator  design  reduces  to  one  of  choosing  the  best  coefficients 
and  then  fabricating  a  filter  which  can  implement  them.  This  chapter  is  a  refinement 
of  work  I  first  reported  in  [65]  on  designing  FIR  filters  for  PM/FM-DD  links.  Links  are 
implemented  using  different  discriminator  Liters,  and  their  performance  is  analyzed 
using  a  small  signal  model,  a  full  signal  model,  and  a  numerical  simulation. 

3.1  Filter  coefficients 

Finite  impulse  response  (FIR)  filters,  with  all  zeros  and  no  poles  in  their  z-transform 
representations,  may  work  well  as  FM  discriminators  because  symmetric  FIRs  can  be 
designed  to  have  exactly  linear  phase,  and  the  theory  shows  that  the  filter’s  phase- 
linearity  affects  the  link’s  linearity.  In  this  and  following  sections,  I  present  sets  of 
FIR  coefficients,  chosen  with  different  criteria,  and  compare  their  performance  as 
discriminators  in  photonic  links. 

My  initial  comparison  is  made  between  different  10th  order  (or  length  11)  sym¬ 
metric  FIR  filters.  The  transfer  function  for  the  positive  slope  filter  goes  from  0  to  1 
within  half  the  Liter’s  free  spectral  range  (FSR),  which  is  the  domain  of  normalized 
angular  frequencies  from  0  to  n.  The  transfer  function  for  the  complemantary  Liter 
with  negative  slope  goes  from  1  to  0  over  the  same  domain.  The  optical  carrier  is 
biased  at  the  midband  angular  frequency  n/2,  which  is  half-held  bias. 

I  chose  three  sets  of  Liter  coefficients.  The  Lrst  two  were  chosen  using  an  optimiza¬ 
tion  routine  with  least-squares  error  minimization.  Because  it  is  difficult  to  match 
the  transfer  function  over  the  full  range,  the  Lrst  Liter  was  optimized  from  0.3  to  0.7. 
The  second  Liter  was  optimized  closer  to  the  carrier  from  0.45  to  0.55.  The  third  set 
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Table  3.1:  Filter  coefficients  for  negative  slope  and  positive  slope,  midband  optimized, 
10th  order,  FIR  discriminators.  Each  filter  is  symmetric,  so  half  the  coefficients  are 
duplicated.  The  symmetric  filters  are  guaranteed  to  have  linear  phase.  The  first  least 
squares  fit  is  optimized  for  normalized  frequencies  0.3  to  0.7,  and  the  second  least 
squares  fit  is  optimized  for  normalized  frequencies  0.45  to  0.55.  The  coefficients  for 
the  maximally  linear  filter  are  from  the  cited  reference.  All  three  filters  are  Type  I 
linear  phase  FIR  filters  (odd-length,  symmetric). 
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of  filter  coefficients  was  chosen  using  the  maximally  linear  criterium.  This  criterium 
was  developed  by  B.  Kumar  and  S.C.  Dutta  Roy  in  [66-68]  for  application  in  digital 
differentiator  filters.  The  maximally  linear  criterium  fixes  a  number  of  derivatives 
of  the  transfer  function  at  a  chosen  frequency,  guaranteeing  high  accuracy  around  a 
small  frequency  band.  If  this  band  is  comparable  to  the  bandwidth  of  modulation, 
overall  I  expect  high  linearity.  The  intuition  for  these  choices  were  based  on  the  error 
bounds  in  the  derived  masks,  which  has  tighter  constraints  close  to  the  carrier. 

The  three  sets  of  filter  coefficients  are  presented  in  Table  3.1  on  page  42.  The 
transfer  functions  for  the  filters  are  plotted  in  Figure  3.1  on  page  43.  All  three  filter 
designs  appear  very  linear  on  the  full  scale,  except  for  the  curvature  at  the  frequencies 
furthest  away  from  the  carrier.  The  figure  also  shows  the  deviation  of  the  transfer 
functions  from  the  ideal  linear  ramp  plotted  on  a  logarithmic  scale.  For  reference, 
I  show  the  cubic  curvature  masks  for  £^max  =  0.01,  0.001,  and  0.0001.  The  first 
least-squares  fit  is  optimized  over  a  wider  range  of  frequencies,  but  the  second-fit 
has  much  smaller  deviation  closer  to  the  carrier.  The  maximally-linear  fit  has  the 
smallest  bandwidth  that  is  optimized,  but  it  is  the  closest  to  the  ideal  filter  over 
that  bandwidth.  This  observation  suggests  a  tradeoff  of  linearity  and  bandwidth  so 
that  the  filter  coefficients  can  be  adjusted  to  best  serve  the  modulation  frequencies  of 
interest. 

The  transfer  functions  were  analyzed  using  the  small  signal  model  presented  in 
the  previous  chapter.  The  code  is  included  in  the  appendix,  Section  A.l.  Fig.  3.2 
shows  OIP3  versus  modulation  frequency  for  links  using  each  of  the  three  sets  of 
filters.  The  link  is  more  linear  for  lower  modulation  frequencies,  and  gets  worse  for 
large  modulation  frequencies.  The  OIP3  from  the  maximally  linear  filter  is  smoothly 
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Figure  3.1:  Transfer  functions  for  the  FIR  discriminators  optimized  at  midband. 
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Figure  3.2:  Simulated  0IP3  for  the  three  different  10th  order  FIR  filter  sets  optimized 
at  midband  versus  normalized  modulation  frequency.  The  photocurrent  is  scaled  for 
10  mA  total  photocurrent  (5  mA  per  detector).  The  filter  is  more  linear  for  lower 
modulation  frequencies,  and  gets  worse  for  large  modulation  frequencies.  For  the  least 
squares  fit  filters,  the  local  minima  for  certain  modulation  frequencies  are  apparent 
in  the  plot. 


Normalized  modulation  frequency 


varying  and  monotonically  decreasing  as  the  modulation  frequency  increases.  These 
are  desirable  properties  when  doing  further  trade-off  analysis  for  the  discriminators, 
so  later  sections  will  employ  the  maximally-linear  filters. 

In  these  designs,  there  is  no  second-order  nonlinearity  for  the  link.  However,  if 
the  detection  is  not  perfectly  balanced,  then  there  is  some  second-order  distortion.  In 
the  simulation,  I  find  that  OIP2  does  not  depend  on  the  modulation  frequency,  but 
does  scale  with  the  common-mode  rejection  ratio  (CMRR)  of  the  detection.  Fig.  3.3 
shows  the  OIP2  versus  CMRR. 


3.2  Scaling  with  filter  order 

Next,  I  study  the  scaling  with  filter  order  of  the  linearity  of  a  link  employing  maximally- 
linear  filters  as  the  discriminator  filters.  Table  3.2  on  page  45  gives  the  coefficients  for 
2nd,  6th,  10th,  14th,  and  18th  order  maximally  linear  filters  optimized  at  half-band 
[67].  The  distortion  of  a  link  using  each  of  the  filters  was  simulated.  Fig.  3.4  shows 
the  scaling  of  OIP3  with  the  filter  order  versus  modulation  frequency.  As  expected, 
higher  order  filters  give  larger  OIP3  than  lower-order  filters.  The  linearity  has  a  large 
improvement  over  a  balanced  MZI  for  a  given  photocurrent.  In  a  physical  link,  the 
highest  values  of  OIP3  (>40  dBm)  may  be  limited  by  photodetector  nonlinearities. 
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Figure  3.3:  Simulated  0IP2  for  the  10th  order  maximally  linear  FIR  filter  set  op¬ 
timized  at  midband  versus  common  mode  rejection  ratio.  The  CMRR  is  given  in 
decibels  of  current  suppressed.  The  photocurrent  is  scaled  for  10  mA  total  photocur¬ 
rent  (5  mA  per  detector).  The  normalized  modulation  frequency  is  0.03,  but  no 
dependence  of  OIP2  on  modulation  frequency  was  observed.  For  infinite  CMRR,  the 
OIP2  value  was  limited  by  the  numerical  precision  of  the  calculation. 


Table  3.2:  Filter  coefficients  for  the  2nd,  6th,  10th,  14th,  and  18th  order  maximally 
linear  filters  in  z-transform  representation.  Each  filter  is  symmetric,  so  half  the  co¬ 
efficients  are  duplicated.  The  coefficients  given  are  for  the  positive  slope  filters.  For 
negative  slope  filters,  the  even-numbered  coefficients  have  opposite  sign. 
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Figure  3.4:  Simulated  0IP3  for  maximally  linear  FIR  filters,  of  different  order,  op¬ 
timized  at  midband  versus  normalized  modulation  frequency.  The  photocurrent  is 
scaled  for  10  mA  total  photocurrent  (5  mA  per  detector). 
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To  study  the  scaling  of  the  spurious-free  dynamic  range  with  filter  order,  I  give 
example  parameters  for  a  phase-modulated  link  based  on  currently  available  commer- 
ical  components.  The  full  free-spectral  range  of  the  filter  is  chosen  to  be  200  GHz, 
and  the  modulation  frequency,  5  GHz.  (The  normalized  frequency  is  therefore  0.05). 
The  transmitter  consists  of  an  external-cavity  laser  with  a  100  kHz  linewidth,  and  a 
lithium  niobate  phase-modulator  with  halfwave  voltage  —  3  V.  The  input  and 
output  loads  are  assumed  to  be  50  fh  The  balanced  detector  handles  up  to  50  mA  of 
dc  photocurrent,  or  25  mA  per  detector. 

Simulating  these  values,  we  end  up  with  a  link  with  gain  of  — 5.6dB  and  noise 
figure  of  10.5  dB.  The  gain  and  noise  figure  will  be  similar  for  a  link  with  a  MZI  dis¬ 
criminator.  The  gain  and  noise  figure  can  be  improved  by  increasing  the  modulation 
efficiency  (decreasing  Vn),  or  increasing  the  detector  power  handling.  The  resulting 
SFDR  versus  filter  order  is  shown  in  Fig.  3.5.  The  phase  noise  limited  and  shot  noise 
limited  approximations  are  shown  in  the  plot  along  with  the  full  noise  figure  model. 
The  SFDR  is  close  to  being  limited  by  the  laser  phase  noise:  there  is  a  difference  of 
about  2  dB.  The  noise  figure  is  primarily  limited  by  the  rf  link- loss,  and  secondly  by 
the  linewidth  of  the  laser.  Both  would  need  to  be  improved  to  get  to  shot  noise  lim¬ 
ited  performance. We  also  find  that  to  some  extent,  one  can  trade-off  between  SFDR 
and  noise  figure  by  adjusting  the  filter’s  FSR.  This  is  shown  in  Fig.  3.6.  With  a 
fixed  modulation  frequency,  making  the  FSR  larger  reduces  the  gain,  thus  making 
the  NF  worse.  However,  the  filter  is  more  linear  closer  to  the  carrier,  so  the  SFDR 
can  be  improved.  Like  the  MZI,  the  SFDR  for  the  2nd  order  order  filter  is  nearly 
independent  of  its  FSR. 
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Table  3.3:  Simulation  parameters 


Parameter 

Value 

fi  (GHz) 

5.0 

h  (GHz) 

5.005 

Input/Output  Impedances  (fl) 

50 

Modulation  efficiency  ( Vn ) 

3 

Laser  linewidth  (kHz) 

100 

Filter  free  spectral  range  (GHz) 

200 

Optical  power  before  filters  (mW) 

250 

Optical  power  incident  upon  each  detector  (mW) 

31.25 

Photodetector  responsivity  (A/W) 

0.8 

DC  photocurrent  per  detector  (mA) 

25 

Third-order  IMD  (GHz) 

4.995 

Second  harmonic  (GHz) 

10 

The  SFDR  with  the  10th  order  filter,  and  200  GHz  FSR  has  a  SFDR  of  129  dB  in  1 
Hz  bandwidth,  which  is  better  than  the  state-of-the-art  links  appearing  in  literature. 
The  SFDR  increases  by  8  dB  for  every  increase  of  4  for  the  filter  order.  This  suggests 
a  great  benefit  from  photonic  integration:  the  link’s  SFDR  scales  with  the  square  of 
the  filter  order! 


3.3  Numerical  link  simulation 

Finally,  I  compare  the  small-signal  model  to  the  full-signal  model  and  numerical 
link  simulations  to  understand  the  limitations  of  the  small-signal  model.  The  full- 
signal  model  is  the  infinite  summation  given  in  (2.9).  For  the  numerical  model,  the 
signal  at  the  output  of  the  link  is  simulated  by  creating  a  time  domain  waveform, 
e(t)  oc  exp[2Tffct  +  27rq  f  performing  a  fast  Fourier  transfer  (FFT),  weighting 
the  frequency  domain  waveform  by  a  given  filter  transfer  function,  performing  an 
inverse  FFT  and  squaring  the  time  domain  waveform  to  obtain  the  photocurrent. 
The  simulation  process  is  illustrated  in  figure  3.7.  The  code  for  all  three  simulations 
are  included  in  the  Appendix.  The  simulation  code  includes  the  effect  of  imperfect 
common-mode  rejection  from  the  balanced  photodetection. 

The  plot  in  Fig.  3.8  shows  the  link  response  of  a  5  GHz  PM-DD  link  using  the 
link  parameters  discussed  in  the  previous  section.  The  fundamental  and  third-order 
intermodulation  distortion  powers  are  plotted  versus  input  power.  The  noise  floor  is 
calculated  using  the  small-signal  model.  As  can  be  seen  in  the  plot,  the  large-signal 
analytical  model  and  the  numerical  simulation  using  FFTs  closely  track  each-other. We 
find  that  the  numerical  simulation  is  much  more  computationally  efficient,  taking  an 
order  of  magnitude  less  time  to  execute. 

For  small  modulation  power,  the  three  models  match  up.  For  large  modulation 
power,  the  distortion  of  the  link  increases  much  faster  than  the  30  dB  per  decade 
suggested  by  the  small-signal  model.  This  can  be  explained  by  observing  that  in 
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Figure  3.5:  Spurious  free  dynamic  range  versus  filter  order  for  5  GHz  PM-DD  links 
using  maximally  linear  filters  and  200  GHz  FSR.  The  link  parameters  are  given  in 
Table  3.3  on  page  47. 


Figure  3.6:  Spurious  free  dynamic  range  for  5  GHz  PM-DD  links  using  maximally 
linear  filters  for  various  FSR. 

Filter  FSR  (GHz) 
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Figure  3.7:  Numerical  model  of  a  PM-DD  or  FM-DD  photonic  link  with  two  discrim¬ 
inator  filters  and  balanced  detection 


Input  signal: 
i(t)  or  V(t) 


frequency  or  phase  modulation,  the  frequency  deviation  of  the  carrier  increases  with 
modulation  depth,  so  more  optical  power  is  spread  into  higher  order  sidebands.  For 
high  modulation  depths,  most  of  the  optical  power  lies  outside  of  the  range  of  frequen¬ 
cies  for  which  the  filter  is  optimized,  creating  more  distortion  than  for  low  modulation 
depths.  For  the  link  under  discussion,  the  models  for  IMD3  begin  to  deviate  around 
1  mW  of  input  power,  or  more  precisely,  the  two  values  are  3  dB  off  when  the  input 
power  is  1.6  dBm.  This  gives  a  phase  modulation  depth  of  0.4  or  a  frequency  modula¬ 
tion  depth  of  2  GHz.  This  makes  sense  in  the  context  of  Carson’s  bandwidth  rule  for 
frequency  modulation:  the  bandwidth  occupied  by  the  modulated  signal  starts  to  have 
a  noticeable  increase  once  the  modulation  frequency  and  the  frequency  modulation 
depth  are  of  the  same  order  of  magnitude. 

The  faster  than  30  dB-per-decade  increase  in  distortion  power  has  consequences 
for  the  actual  spurious  free  dynamic  range  seen  by  the  system.  Fig.  3.10  compares 
the  small-signal  approximation  to  the  SFDR  with  SFDR  values  calculated  by  finding 
the  intercept  of  the  IMD3  with  the  noise  in  a  given  bandwidth.  Large  bandwidths 
will  be  unequally  affected  by  the  sideband  spill-over  effect.  A  link  designer  needs  to 
be  cognizant  of  the  full-signal  model  in  order  to  accurately  predict  the  spurious  signal 
levels  seen. 


3.4  Summary 

We  have  demonstrated  by  simulation  that  frequency  and  phase  modulated  microwave- 
photonic  links  with  very  high  linearity  are  obtainable  by  using  FIR  optical  filters  to 
perform  demodulation.  Links  using  filters  designed  with  the  maximally  linear  criteria 
greatly  exceed  the  linearity  performance  of  an  MZI.  Although  this  chapter  did  not 
simulate  the  performance  of  HR  filters  as  discriminators,  in  general  it  is  expected 
that  more  closely  matched  filters  could  be  implemented  with  fewer  stages  with  an  HR 
architecture.  We  have  observed  that  linearity  degrades  for  high  modulation  depths  as 
power  is  spread  into  high-order  optical  sidebands  far  from  the  optical  carrier.  We  find 
that  a  tenth-order  FIR  filter  designed  using  the  maximally  linear  criteria  can  obtain  a 
129  dB  •  Hz2/3  SFDR  with  50  mA  of  photocurrent.  The  SFDR  scales  with  the  square 
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Figure  3.8:  Link  response  versus  input  power  for  a  5  GHz  PM-DD  link  using  tenth- 
order  maximally  linear  filters.  The  link  parameters  are  given  in  the  text. 
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Figure  3.9:  Link  response  versus  input  power  for  a  5  GHz  PM-DD  link  using  maxi¬ 
mally  linear  filters  of  different  order. 
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Figure  3.10:  Spurious  free  dynamic  range  versus  bandwidth  for  5  GHz  PM-DD  links 
using  maximally  linear  filters  of  different  orders. 


of  the  FIR  filter  order,  suggesting  a  benefit  to  increasing  photonic  integration. 
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Chapter  4 

Phase  modulation  experiments 


Using  complementary  linear-field  discriminator  filters,  we  believe  we  have  demon¬ 
strated  PM-DD  and  FM-DD  links  with  the  highest  linearities  which  have  been  pub¬ 
lished  thus  far,  as  measured  by  third-order  and  second-order  output-intercept  points 
(OIP3  and  OIP2)  normalized  to  a  fixed,  photodetector-limited  photocurrent.  Our  dis¬ 
criminator  filters  are  fabricated  in  a  low-loss  silica-on-silicon,  planar-lightwave-circuits 
(PLC)  process  at  Alcatel-Lucent  Bell  Laboratories.  We  report  link  measurements  us¬ 
ing  both  a  cascaded  MZI  FIR  lattice  filter  and  a  ring  assisted  MZI  (RAMZI)  HR 
filter,  and  with  both  phase  modulation  and  frequency  modulation. 

The  discriminators  are  based  on  two  architectures:  a  cascaded  MZI  FIR  lattice 
filter  [69]  and  a  ring  assisted  MZI  (RAMZI)  HR  filter  [70].  For  both  types  of  discrimi¬ 
nators,  we  demonstrate  >  6  dB  improvement  in  the  link’s  third-order  output  intercept 
point  (OIP3)  over  a  MZM  link.  We  show  that  the  links  have  low  second-order  dis¬ 
tortion  when  using  balanced  detection.  Using  high  optical  power,  we  demonstrate  an 
OIP3  of  39.2  dBm.  We  also  demonstrate  4.3  dB  improvement  in  signal  compression. 

4.1  Planar  lightwave  circuit  filters 

A  discriminator  filter  approximating  the  ideal  complementary  linear-field  response 
can  be  constructed  using  silica-on-silicon  planar  lightwave  circuits  (PLC)  [71].  The 
transform  function  of  an  FIR  filter  can  be  realized  in  PLC  with  just  MZIs  and  di¬ 
rectional  couplers.  One  implementation  of  a  multi-stage  optical  FIR  filter  in  PLC 
is  the  lattice  filter  [64],  The  lattice  filter  architecture  has  a  low-loss  passband  and 
requires  only  N+l  couplers  for  an  Nth  order  filter,  which  are  advantages  over  other 
optical  filter  architectures.  The  lattice  filter  architecture  is  shown  in  Figure  4.1  on 
page  53,  indicating  for  each  stage  the  coupling  coefficients,  designated  by  k,  and  the 
phase  shifts,  designated  by  ip.  Each  stage  has  a  unit  delay,  z*1 .  The  dashed  lines 
indicate  additional  filter  stages  omitted  from  the  figure.  A  recursion  relation  exists 
that  transforms  between  given  filter  coefficients  and  the  corresponding  coupling  ratios 
and  phase  shifts  [64],  The  recursion  relation  for  the  tenth-order  lattice  filter  design 
gives  210,  or  1024  solutions.  For  a  tenth-order  maximally  linear  discriminator  filter, 
whose  coefficients  were  given  in  the  previous  chapter  in  Table  3.2  on  page  45,  one 
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Figure  4.1:  FIR  lattice  filter  architecture 


(pi  cp2  <plO 


Table  4.1:  Filter  phase  and  coupler  parameters  for  a  tenth-order  maximally  linear 
discriminator  filter  in  lattice  filter  form _ 


Phase  shift 

Value 

Coupling  ratio 

Value 

Tunable  coupler  phase 

Ai 

0 

«  0 

0.674271 

0.607391 

<^2 

—  7T 

«1 

0.635460 

0.648224 

</?3 

0 

0.837472 

0.414954 

<P4 

—  7 T 

0.514751 

0.770645 

V 5 

0 

K4 

0.918512 

0.289487 

<^6 

0 

hi5 

0.538578 

0.746782 

<P7 

—  7T 

0.918384 

0.289721 

Vs 

0 

k7 

0.515217 

0.770179 

<P9 

—7 r 

0.837319 

0.415161 

—7 r 

hig 

0.635828 

0.647842 

K 10 

0.325694 

0.963443 

particular  solution  for  the  parameters  in  lattice  filter  form  is  listed  in  Table  4.1. 

Up  to  tenth-order  FIR  lattice  filters  have  been  implemented  in  PLC  for  various 
applications.  A  research  group  at  NTT  laboratories  has  extensively  explored  tunable 
optical  FIR  lattice  filters.  Tunable  coupling  ratios  are  implemented  by  using  sym¬ 
metric  Mach-Zehnder  interferometers  with  thermal-optic  phase  shifters.  A  diagram 
of  a  tunable  FIR  filter  is  shown  in  Figure  4.2  on  page  53.  The  intended  application 
is  dispersion  compensation,  but  because  the  filters  are  tunable,  they  can  be  used  for 
any  filter  transfer  function  desired,  including  discriminators. 

The  group  has  fabricated  eight-order  filters  in  silica  with  chromium  heaters,  with 
50  GHz  FSR  [69]  and  200  GHz  FSR  [72],  and  arrays  of  fifth-order  filters  with  50 


Figure  4.2:  Tunable  PLC  FIR  lattice  filter  architecture 
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GHz  FSR  [73-75].  They  claim  control  of  the  individual  phase  shifters  to  accuracy 
better  than  0.017T  radians.  For  the  fifth-order  filters,  to  reduce  the  required  bias 
power  on  the  heaters,  they  use  a  phase-trimming  technique  that  involves  introducing 
heating  induced  stress.  They  have  also  proposed  a  100  GHz  FSR  filter  in  a  reflection 
configuration  to  double  its  effective  length  [76]. 

A  collaboration  between  Siemens,  University  of  Kiel  and  IBM  Research  Zurich  has 
implemented  the  same  architecture  on  a  more  compact  silicon  oxynitride  platform  [77]. 
The  applications  include  both  EDFA  gain  equalization  and  dispersion  compensation. 
They  have  demonstrated  sixth,  seventh  and  tenth-order  filters  with  100  GHz  FSR. 
The  collaboration  has  explored  a  number  of  adaptive  feedback  approaches  for  setting 
the  filter’s  phase  shifters  [78]. 

•  Optical  spectrum  analysis  [79]:  They  have  used  an  optical  spectrum  analyzer 
to  compare  the  amplified  spontaneous  emission  spectrum  to  a  desired  intensity 
profile.  A  computer  running  the  Levenberg-Marquart  optimization  algorithm  (a 
modified  Gauss-Newton  algorithm)  varies  the  power  to  the  phase  shifters  until 
the  desired  profile  is  obtained. 

•  Electrical  spectrum  monitoring  [80-82]:  ESM  is  another  feedback  approach, 
where  power  at  certain  frequencies  are  used  as  a  feedback  mechanism.  Pilot 
tones  or  knowledge  about  the  signaling  over  the  link  determines  the  optimal 
choice  of  electrical  filters. 

•  Eye  opening  [83-85]:  An  adaptive  feedback  approach  for  digital  signals  looks  at 
an  eye  diagram  and  uses  a  Levenberg-Marquart  optimization  to  maximize  the 
eye  opening. 

•  LMSE  /  minimize  ISI  [80,  82,  86-88]:  Another  method  for  digital  signals  uses 
minimization  between  the  decision  and  signal  as  a  feedback  signal.  These  meth¬ 
ods  are  not  suitable  for  analog  links. 

For  simplicity  and  cost,  setting  the  filters  coefficients  without  using  a  feedback  system 
is  desired.  The  IBM  collaboration  has  developed  a  calibration  procedure  to  produce 
a  table  look-up  for  tunable  coupler  and  phase  shifter  responses  versus  applied  tuning 
power  [79].  The  technique  uses  the  OSA  approach  to  iteratively  tune  all  couplers  to 
zero  cross  coupling.  There  is  a  procedure  to  individually  characterize  each  tunable 
coupler  and  asymmetric  MZI  by  measuring  the  output  power  versus  tuning.  The 
filter  then  can  be  set  to  a  pre-calculated  inverse  system.  An  alternative  calibration 
approach  is  given  by  the  NTT  group  in  [89]  that  does  not  require  a  feedback  loop. 
The  approach  uses  incoherent  light  to  characterize  each  tunable  coupler,  and  coherent 
light  to  characterize  the  asymmetric  Mach  Zehnders. 

4.2  Implementation  and  characterization 

As  part  of  this  program  of  research,  two  types  of  filters  were  fabricated  and  packaged 
at  Alcatel-Lucent  Bell  Laboratories  by  Dr.  Mahmoud  Rasras:  a  cascaded  MZI  FIR 


54 


Figure  4.3:  (a)  Filter  stage  for  an  FIR  lattice  filter  (b)  Filter  stage  for  an  HR,  RAMZI 
filter. 


lattice  filter  and  a  ring  assisted  MZI  (RAMZI)  IIR  filter.  A  single  stage  of  each 
filter  is  illustrated  in  Fig.  4.3.  These  filters  can  be  thermally  tuned  using  chromium 
heaters  to  implement  arbitrary  filter  transfer  functions.  The  RAMZI  IIR  filter  is  a 
third-order  filter  with  an  all  pass  ring  resonator  structure  coupled  to  the  delay  arm 
of  an  MZI.  The  FIR  filter  is  a  sixth-order  filter  with  120  GHz  free-spectral  range. 
Our  filter  has  6  stages  of  symmetrical  MZIs  (switches)  and  asymmetrical  MZIs  (delay 
line  interferometers)  which  are  tunable  using  chromium  heaters  deposited  on  the 
waveguides.  Figs.  4.4  and  4.5  show  photographs  of  a  fabricated  and  packaged  FIR 
filter. 

The  experimental  system  for  a  phase-modulated  link  measurement  is  illustrated  in 
Figure  4.6  on  page  57.  A  polarization  tracker  is  used  at  the  output  of  the  ECTL,  and, 
where  possible,  the  optical  paths  are  polarization  maintaining  fiber.  Two  tunable  RF 
sources  are  combined  to  modulate  a  commercial  lithium  niobate  phase  modulator  to 
perform  two-tone  distortion  measurements.  We  use  a  personal-computer-based  analog 
output  card  to  generate  bias  currents  for  the  heaters  to  tune  the  transfer  function  of 
the  discriminator.  See  Figs.  4.7  and  4.8.  The  paths  between  the  filters  and  balanced 
detectors  are  trimmed  to  match  delay  and  attenuation.  For  the  FM  measurements, 
the  tunable  laser,  polarization  controller,  and  phase  modulator  are  replaced  with  the 
directly  modulated  FM  laser. 

4.3  Link  Results 

4.3.1  Phase-modulated  link  with  FIR  filter 

We  performed  link  measurements  using  the  FIR  filter  and  phase  modulation.  In 
our  experiment  the  discriminator  filter  is  dynamically  tuned  to  minimize  the  link 
distortion.  The  filter  has  13  degrees-of-freedom  to  adjust.  If  the  filter  is  ideal,  one  can 
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Figure  4.4:  Photograph  of  single  FIR  filter  with  wiring  board  inside  protective  box. 


Figure  4.5:  Photograph  of  single  FIR  filter  mounted  on  heat  sink. 
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Figure  4.6:  Diagram  of  the  system  used  for  characterization 
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Figure  4.7:  Photograph  of  current  amplifier  board  to  drive  the  chrome  heaters  on  the 


tunable  filters 
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Figure  4.8:  Photograph  of  National  Instruments  analog  input/output  card  interface. 


in  principle  choose  all  the  parameters  a-priori  to  implement  desired  filter  coefficients. 
However,  it  is  difficult  to  characterize  precisely  the  correspondence  between  currents 
applied  to  each  waveguide  heater  and  the  resulting  optical  phase  shift.  Imperfections 
in  the  filter  fabrication  also  make  the  characterization  difficult.  Therefore,  feedback 
is  used  to  choose  the  correct  biases  to  the  heaters. 

We  use  an  optimization  routine  employing  a  downhill-simplex  algorithm  to  tune 
the  heater  settings  for  the  discriminator  filter.  Two  radio  frequency  synthesizers  are 
used  to  generate  tones  at  2  GHz  and  2.0001  GHz  with  equal  RF  powers.  The  error 
signal  for  the  optimization  routine  is  the  third-order  intermodulation  distortion  term 
at  1.9999  GHz,  normalized  to  the  dc  photocurrent  and  the  fundamental  signal  power. 
The  start  point  for  each  heater  is  randomly  chosen  within  an  acceptable  range  of 
currents  which  will  not  cause  damage  to  the  device.  The  routine  varies  the  heater 
settings  to  minimize  the  error  signal,  thus  maximizing  the  OIP3.  The  routine  reaches 
a  minimum  error  value  after  less  than  100  iterations. 

One  of  the  filters  was  tuned  to  the  desired  linear  ramp  and  linear  phase  transfer 
function.  The  phase  and  amplitude  of  the  filter  were  measured  with  an  optical  vector 
network  analyzer  (OVNA).  The  transfer  function  shown  in  Fig.  4.9  is  normalized 
to  a  7  dB  filter  insertion  loss.  The  insertion  loss  could  be  improved  by  better  fiber 
coupling  into  the  filter.  The  waveguide  loss  for  silica  PLC  is  not  a  significant  loss 
mechanism.  At  the  50%  field  amplitude  transmission  point,  both  the  amplitude  and 
phase  of  the  transfer  function  appear  linear  within  the  accuracy  of  the  instrument. 

We  report  distortion  measurements  made  with  a  single  branch  of  the  filter  and 
single-ended  detection.  With  tones  at  2  GHz  and  2  GHz  +  100  kHz  for  the  funda¬ 
mental  frequencies,  we  stepped  the  wavelength  of  the  laser  to  determine  the  optimal 
bias  point  on  the  filter.  At  each  wavelength,  we  collected  the  receiver  power  at  2  GHz 
and  the  third-order  intermodulation  distortion  (IMD3)  power  at  2  GHz  -  100  kHz. 
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Figure  4.9:  Achieved  filter  amplitude  and  phase  for  the  6th  order  FIR  lattice  filter. 


Wavelength  (nm) 

Fig.  4.10  shows  the  fundamental  power  and  third-order  nonlinearity  as  a  function 
of  carrier  wavelength.  The  distortion  remains  low  over  a  wavelength  span  of  10  pm. 
A  1549.937  nm  carrier  wavelength  gives  an  optimal  ratio  of  fundamental  to  IMD3 
power.  The  optimal  operating  wavelength  corresponded  with  a  point  where  the  filter 
has  50%  amplitude  transmission.  For  a  system  with  a  fixed  wavelength  of  operation, 
the  filter  itself  can  be  tuned  in  wavelength  by  adjusting  the  phase  delays  in  each 
stage.  A  single  filter  could  also  be  tuned  to  accommodate  a  variety  of  sources  at 
different  wavelengths,  such  as  wavelength  division  multiplexing  channels,  since  the 
filter  transfer  function  repeats  over  each  free-spectral  range. 

At  the  optimal  wavelength,  we  varied  the  RF  power  input  into  the  link  and  mea¬ 
sured  the  IMD3  and  fundamental  power.  The  data  is  shown  in  Fig.  4.11.  The 
distortion  clearly  showed  a  cubic  dependence  with  input  power.  For  a  photocurrent 
of  0.11  mA,  we  measured  an  OIP3  of  -19.5  dBm.  The  OIP3  of  a  PM  link  using  an 
MZI  discriminator  and  the  same  photocurrent  is  -26.2  dBm.  This  particular  transfer 
function  displayed  a  6.7  dB  OIP3  performance  improvement  over  an  MZI  with  the 
sained  received  photocurrent.  For  shot-noise  limited  noise  performance,  the  link  has 
a  spurious  free  dynamic  range  of  112  dB  ■  Hz2//3.  If  the  light  is  amplified  to  produce  10 
mA  of  photocurrent,  OIP3  increases  to  19.7  dBm  and  the  shot-noise  limited  spurious 
free  dynamic  range  is  125  dB  •  Hz2/3. 

These  initial  results  using  optical  lattice  filters  to  discriminate  phase  modulation 
show  appreciable  improvement  in  third-order  nonlinearity  over  the  MZI,  achieving  6.7 
dB  improvement  in  OIP3  at  2  GHz  modulation  frequency,  for  a  10  pm  (1.25  GHz) 
span  of  carrier  wavelengths.  Because  electrical  spectrum  monitoring  and  a  feedback 
algorithm  are  used  to  optimize  the  filter  coefficients,  the  system  is  robust  to  fabri¬ 
cation  imperfections  and  other  variations  in  the  devices.  Our  particular  FIR  lattice 
filter  had  variations  in  switch  contrast  and  waveguide  loss  that  were  corrected  by  the 
optimization.  No  assumptions  were  made  about  the  mapping  of  the  heater  settings 
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Figure  4.10:  Fundamental  and  third-order  intermodulation  distortion  versus  laser 
wavelength.  The  modulation  power  is  fixed  at  10  dBm  and  the  photocurrent  is  fixed 
at  0.11  mA. 


Figure  4.11:  Fundamental  and  third-order  intermodulation  distortion  versus  modu¬ 
lation  power.  The  photocurrent  is  fixed  at  0.11  mA  and  the  wavelength  is  fixed  at 
1593.7  nm. 
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Figure  4.12:  Achieved  filter  amplitude  and  phase  for  the  RAMZI  filter. 


Wavelength  (nm) 


for  the  filter  to  its  transfer  function  to  run  the  optimization  routine,  so  the  technique 
could  be  used  for  any  optical  filter  architecture.  Theory  and  simulation  suggest  that 
much  larger  improvements  in  linearity  are  possible  through  careful  adjustment  of  the 
filter  transfer  function. 

Besides  just  further  improvement  in  nonlinearity,  the  optimization  goal  can  also  be 
changed  to  seek  improvement  in  signal  gain,  or  nonlinearity  improvement  over  chosen 
carrier  wavelengths  or  modulation  frequency  spans.  High-order  filters  could  also  be 
optimized  to  correct  for  gain  equalization  and  fiber  dispersion.  For  these  reasons, 
using  general  purpose,  tunable  filters  for  phase  discrimination  is  a  very  promising 
technique  and  suitable  for  much  further  study. 

4.3.2  Phase-modulated  link  with  HR  filter 

We  also  performed  PM-DD  measurements  with  the  HR  filter.  Fig.  4.12  shows  the 
tuned  transfer  function  of  one  of  the  branches  of  the  RAMZI  filter.  The  transfer 
function  was  tuned  by  hand  using  the  OVNA  measurement  for  reference.  The  plotted 
transfer  function  is  normalized  to  5  dB  insertion  loss.  The  transfer  function  was 
measured  with  an  OVNA.  The  second  branch  was  tuned  to  a  transfer  function  with 
opposite  slope. 

For  our  first  distortion  measurement,  we  use  fundamental  frequencies  of  5  GHz  and 
5  GHz  +  100  kHz.  Like  with  the  FIR  filter,  we  stepped  the  wavelength  of  the  tunable 
laser  and  measured  the  fundamental  and  IMD3,  but  in  this  case,  the  measurement  is 
performed  with  the  balanced  detection  and  two  Liters.  The  optical  power  is  amplified 
so  that  the  total  DC  photocurrent  added  from  the  two  detectors  is  10.5  mA..  Fig. 
4.13  shows  the  calculated  OIP3  from  the  fundamental  and  IMD3  data.  Also  shown 
on  the  graph  is  the  theoretical  OIP3  from  a  dual-output  MZM  with  the  same  DC 
photocurrent,  which  we  exceed  for  a  range  of  biases.  It  is  important  to  note  that 
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Figure  4.13:  Output  intercept  point  of  third-order  intermodulation  distortion  versus 
laser  wavelength  in  simulation  and  experiment.  The  total  photocurrent  is  fixed  at 
10.5  mA  and  the  modulation  frequency  is  5  GHz.  The  theoretical  OIP3  of  a  link  with 
a  dual-output  MZM  and  the  same  received  photocurrent  is  also  plotted  in  the  figure. 


Wavelength  (nm) 


the  measured  and  theoretical  power  values  include  a  6  dB  power  loss  due  impedance 
matching  for  the  photodetector 

Furthermore,  the  measured  OVNA  transfer  function  was  used  in  a  numerical  link 
simulation,  as  described  in  Chapter  3.  The  simulation  trend  of  OIP3  versus  wave¬ 
length  matched  with  the  experimental  data,  except  for  peaks  around  the  optimal  bias 
points.  This  may  be  due  to  optical  system  factors  not  included  in  the  simulation, 
such  as  balancing  the  two  filters,  fiber  dispersion,  polarization  drift,  back  reflections, 
etc.,  or  may  be  due  to  lack  of  resolution  of  the  OVNA  in  measuring  ripples  in  the 
transfer  function. 

We  tuned  the  wavelength  to  1549.964  nm,  which  is  approximately  50%  field  trans¬ 
mission  bias  and  around  the  maximum  measured  OIP3  at  5  GHz.  Varying  the  modula¬ 
tion  frequency,  we  measured  both  the  IMD3  and  the  second  harmonics,  and  calculated 
OIP3  and  OIP2,  shown  in  Fig.  4.14.  For  frequencies  between  1  to  10  GHz,  we  exceed 
the  OIP3  of  the  MZM.  The  best  frequency  was  around  4  GHz,  which  gave  us  an  8.5 
dB  OIP3  improvement  over  the  dual-output  Mach-Zehnder.  Because  it  is  a  balanced 
device,  the  second-order  distortion  is  also  low  over  the  whole  band.  The  link  gain 
increases  for  higher  frequencies  because  of  the  1/f  PM  to  FM  correspondence. 

We  compared  the  PM  link  and  a  dual-output  MZM  link  using  a  high-power  bal¬ 
anced  photodetector  array  in  development  [90].  As  shown  in  Fig.  4.15,  the  PM-DD 
link  achieved  4.3dB  better  RF  compression  power  than  the  dual  MZM.  The  link  is 
nearly  transparent,  with  only  3.8  dB  signal  loss.  As  shown  in  Fig.  4.16,  increasing 
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Figure  4.14:  0IP3  and  0IP2  versus  modulation  frequency  at  a  fixed  photocurrent  of 
10.5  mA  and  wavelength  of  1549.964  nm. 


RF  Frequency  (GHz) 

the  effective  photocurrent,  we  achieved  a  maximum  OIP3  of  39.2  dBm  at  4  GHz 
modulation  frequency. 

4.4  Summary 

In  this  chapter,  I  have  reviewed  previous  work  in  implementing  FIR  filters  in  photonics 
using  optical  lattice  filters,  and  discussed  experimental  results  using  phase  modula¬ 
tion.  Using  a  sixth  order  optical  lattice  filter  and  a  feedback  optimization  system,  we 
instantiated  a  linear-field  and  linear-phase  filter  transfer  function.  This  gave  a  6.7  dB 
improved  OIP3  over  an  MZI  for  a  2  GHz  phase  modulated  link.  Using  a  3rd-order 
HR  filter,  we  demonstrated  improved  OIP3  over  an  MZI  for  frequencies  between  1  to 
10  GHz,  with  the  optimal  improvement  of  8.5  dB  at  frequency  4  GHz.  Using  high 
optical  power,  we  demonstrate  an  OIP3  of  39.2  dBm.  The  link  had  low  second-order 
distortion  when  using  balanced  detection. 
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Figure  4.15:  Output  power  versus  modulation  power  compared  to  a  dual-output 
Mach-Zehdner  modulator  measured  experimentally.  The  frequency  is  fixed  at  3.3 
GHz  and  the  effective  DC  photocurrent  at  141  m A. 


Figure  4.16:  OIP3  versus  effective  DC  photocurrent.  The  frequency  is  fixed  at  4.0 
GHz  and  the  modulation  power  at  0  dBm. 
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Chapter  5 

Frequency  modulation  experiments 


In  this  chapter,  I  present  a  review  of  work  in  the  field  of  FM  lasers,  the  experimental 
characterization  of  current  modulated,  distributed  Bragg  reflector  lasers  developed 
for  this  program  of  research,  and  microwave  photonic  link  results  using  the  FM  lasers 
and  PLC  discriminator  filters. 


5.1  Review  of  FM  lasers 

There  has  been  much  work  over  the  years  in  designing  frequency  modulated  lasers. 
For  application  in  microwave  photonic  links,  the  theory  in  Chapter  2  indicates  that 
the  FM  lasers  must  have  small  linewidths  (<  100  kHz)  and  large  modulation  effi¬ 
ciency  (of  the  order  of  10s  of  GHz/mA  or  100s  of  GHz/V)  in  order  to  give  links 
with  small  noise  figures.  In  addition,  the  lasers  also  must  have  low  residual  intensity 
modulation  (Section  2.6.3)  and  low  intrinsic  distortion  in  order  to  ensure  linearity. 
Frequency  modulation  can  be  produced  in  semiconductor  lasers  by  active  gain  mod¬ 
ulation,  phase  modulation,  or  grating  modulation  of  single  or  multi-section  devices. 
Both  distributed  Bragg  reflector  (DBR)  and  distributed  feedback  (DFB)  lasers  have 
been  used  in  frequency  modulation  experiments.  The  phase  shift  in  each  section  is 
either  produced  by  current  or  voltage  modulation.  In  current  modulation,  the  junc¬ 
tion  is  forward  biased,  and  injected  carriers  change  the  refractive  index  through  the 
free-carrier  plasma  effect.  In  voltage  modulation,  the  junction  is  reverse  biased,  and 
a  field  effect  such  as  the  Franz-Keldysh  (FK)  effect  in  bulk  material  or  the  quantum- 
confined  Stark  effect  (QCSE)  in  quantum-well  material  produces  a  refractive  index 
change.  For  a  low  efficiency  approach  to  frequency  modulation,  [91]  and  [92,  93]  have 
suggested  using  phase  modulators  operated  past  their  RC  limited  bandwidth. 

There  have  been  a  number  of  theoretical  papers  looking  at  laser  distortion  under 
frequency  modulation.  [94-96]  used  Bessel  function  and  perturbation  models  to  solve 
for  harmonic  and  intermodulation  distortion  in  laser  IM  and  FM.  In  [48,  97-102],  the 
authors  used  Taylor-series  analyses  to  estimate  the  harmonic  and  intermodulation 
distortion  of  an  FM  laser  after  a  Mach  Zehnder  interferometer.  Unfortunately,  none 
of  these  papers  applied  their  theory  to  multi-section  devices. 


65 


5.1.1  Fabry-Perot  lasers 

In  some  of  the  earliest  work  on  frequency-modulation  of  semiconductor  lasers,  [42] 
measured  the  direct  FM  and  IM  characteristics  of  three  single-section  Fabry-Perot 
AlGaAs  lasers.  The  results  were  consistent  with  models  for  thermal  and  carrier  injec¬ 
tion  effects.  The  device  response  not  flat  and  decreased  several  orders  of  magnitude 
from  low  to  high  modulation  frequency.  Modulation  efficiencies  were  of  the  order  of  1 
GHz/mA  for  low  frequencies  (thermal),  and  100  MHz/mA  for  high  frequencies  (car¬ 
rier  injection),  which  is  not  sufficient  for  microwave  photonic  links.  The  nonlinearity 
of  the  modulation  was  not  studied. 

Later  Fabry-Perot  devices  used  a  second  tuning  section  and  QCSE  based  voltage 
modulation  to  increase  the  modulation  bandwidth.  The  device  reported  in  [103] 
had  a  flat  response  up  to  2  GHz,  frequency  modulation  efficiency  of  7  GHz/V,  and 
3.7  MHz  minimum  linewidth.  [104]  claimed  a  simpler  fabrication  technique  for  the 
same  type  of  device.  Their  improved  laser  had  flat,  RC  limited  response  from  5 
kHz  to  500  MHz,  with  a  modulation  efficiency  of  20  GHz/V.  In  addition,  the  laser 
also  had  very  low  residual  IM.  The  residual  intensity  modulation  depth  was  0.05  for  a 
frequency  deviation  of  4  GHz.  To  demonstrate  the  benefit  of  the  reverse  biased,  QCSE 
modulation,  measurements  were  compared  with  forward  biased  carrier  injection.  The 
efficiency  ranged  from  1.7  GHz/mA  at  5  kHz  to  0.3  GHz/mA  at  50  MHz.  In  a 
later  device,  [105],  the  capacitance  was  improved,  thus  giving  a  uniform  FM  response 
from  30  kHz  to  6  GHz,  and  FM  efficiency  of  about  10  GHz/V.  The  linewidth  was 
approximately  10  MHz.  The  residual  IM  was  low  for  this  device:  the  IM  depth  was 
0.05  for  a  frequency  deviation  of  4  GHz.  For  this  device  under  forward  bias,  the  FM 
efficiency  was  1  GHz/mA  at  low  frequency  but  fell  off  rapicly  to  10  MHz/mA  at  1 
GHz.  The  residual  IM  depth  was  0.09  for  frequency  deviation  of  3.8  GHz. 

5.1.2  DBR  lasers 

Multi-section  distributed  Bragg  reflector  lasers  have  also  been  explored.  In  [106],  the 
authors  presented  a  two  section  DBR  laser  with  grating  section  modulation  and  cur¬ 
rent  injection.  The  modulation  bandwidth  was  400  MHz  and  the  modulation  efficiency 
was  3  GHz/mA.  The  linewidth  was  approximately  50  MHz.  For  1  GHz  frequency  de¬ 
viation,  the  RIM  depth  was  0.01.  Another  author,  [107],  used  a  multi-electrode  DBR 
laser  to  produce  pure  frequency  modulation.  Because  carrier  injection  in  a  passive 
grating  section  has  limited  bandwidth,  they  instead  modulated  the  active  region  of 
the  laser  with  a  non  uniform  current  density.  With  only  single  electrode  modulation, 
the  bandwidth  was  600  MHz  ,  the  modulation  efficiency  was  200  MHz/mA,  and  the 
IM  efficiency  was  0.014/mA.  The  bandwidth  was  limited  by  the  capacitance  of  the  de¬ 
vice.  With  push-pull,  dual  electrode  modulation,  the  FM  efficiency  was  360  MHz/mA 
and  IM  efficiency  was  0.00252/mA.  In  [108],  a  5  MHz  linewidth  was  reported  for  the 
same  device.  [109]  carried  out  a  theoretical  analysis  of  current  injected,  three-section 
(gain,  phase,  grating)  DBR  lasers  using  general  rate  equations.  The  authors  solved 
for  the  frequency  dependence  of  the  FM  and  IM  modulation  efficiencies. 

Later  DBR  devices  used  voltage  modulation.  The  laser  in  [110]  used  a  MQW  gain 
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section  and  superlattice  passive  sections.  The  authors  measured  a  4  GHz  bandwidth, 
with  38.75  GHz/V  modulation  efficiency.  A  similar  device  was  reported  by  [111,  112], 
which  used  the  Franz-Keldysh  effect  to  tune  the  Bragg  section.  The  linewidth  was 
less  than  10  MHz  over  the  tuning  range.  The  authors  published  another  device  using 
an  electro-absorption  modulated  section  in  [113,  114].  [115]  theoretically  studied  how 
to  optimize  phase-section  modulated  QCSE  DBR  lasers.  The  theoretical  modulation 
efficiency  was  2.5  GHz/V,  and  the  spurious  intensity  modulation  was  suppessed.  More 
recently,  [116]  presented  a  three-section  InGaAsP  DBR  laser.  In  this  device,  the 
passive  and  active  use  the  same  MQW  material,  but  the  absorption  was  reduced  in  the 
active  sections  by  using  quantum-well  intermixing  (QWI)  based  on  ion  implantation. 
The  modulation  response  was  uniform  from  10  MHz  to  10  GHz.  The  tuning  efficiency 
of  the  DBR  section  was  375  GHz/V  at  dc.  In  [117,  118],  the  authors  modulated  the 
phase-section  of  a  DBR  laser.  To  reduce  residual  IM  and  increase  the  bandwidth,  the 
authors  designed  the  bandgap  of  the  phase  control  region  to  be  much  wider  than  that 
of  the  gain  region,  so  the  waveguide  loss  variation  induced  by  the  refractive  index 
modulation  would  be  extremely  small.  The  authors  presented  an  array  of  four  lasers 
in  [119].  The  modulation  efficiency  was  approximately  5  GHz/V.  The  linewidth  of 
the  lasers  was  15  MHz. 

5.1.3  DFB  lasers 

Finally,  authors  have  studied  FM  distributed  feedback  lasers  with  single  and  multiple 
electrodes.  [120]  measured  the  modulation  response  versus  frequency  of  a  commerical 

1.3  micron  DFB  laser.  The  bandwidth  was  approximately  4.3  GHz,  but  it  was  not 
flat,  there  was  an  observeable  resonance  peak.  At  1  GHz,  the  FM  response  was  170 
MHz/rnA,  and  the  IM  response  was  0.022  /  mA.  [121]  developed  a  two  section  DFB 
laser,  with  a  DFB  active  region  and  a  phase  tuning  region.  The  phase  region  was 
modulated  with  200  MHz  bandwidth.  The  maximum  FM  efficiency  at  10  MHz  was 
16  GHz/mA.  The  linewidth  was  of  the  order  of  50  MHz.  The  intensity  modulation 
was  less  than  0.01  for  1  GHz  frequency  deviation.  [122]  reported  a  complex  coupled 
DFB  laser  with  a  specially  engineered  grating.  It  had  a  flat  FM  response  from  10  kHz 
to  more  than  20  GHz,  limited  by  the  measurement  system.  The  modulation  efficiency 
was  0.95  GHz/mA  at  1  GHz,  and  the  linewidth  was  12  MHz. 

[29,  123]  gave  a  theoretical  model  for  a  two  section  DFB  laser  using  a  QCSE  phase 
section.  In  the  results,  the  highest  FM  efficiency  reported  was  80  GHz/V  at  about 
3  GHz  modulation.  The  modulation  was  not  flat,  but  the  bandwidth  was  greater 
than  several  GHz.  The  authors  argued  that  DBR  lasers  are  better  than  DFB  lasers 
because  the  phase-modulating  section  can  occupy  a  relatively  larger  fraction  of  the 
cavity  leading  to  higher  FM  efficiency. 

[124]  presented  a  multi-electrode  DFB  device.  Multiple  electrodes  were  engineered 
so  that  carrier  and  thermal  effects  would  be  in  phase.  There  were  electrodes  on  the 
center  and  sides  of  the  active  DFB  region.  The  modulation  efficiency  was  flat  from  100 
kHz  to  more  than  15  GHz.  At  100  kHz,  the  modulation  efficiency  was  0.64  GHz/mA. 
At  1  GHz,  it  was  0.37  GHz/mA.  The  authors  used  a  long  1.2  mm  cavity  to  reduce 
the  linewidth  to  less  than  1  MHz.  [125]  reported  a  DFB  laser  with  5  electrodes.  The 
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Figure  5.1:  Self  heterodyne  laser  linewidth  measurement  experimental  setup. 


bandwidth  was  4  GHz,  the  FM  efficiency  was  0.65  GHz/mA,  and  the  linewidth  was 
3  MHz. 

[30,  126,  127]  reported  a  high  efficiency  FM  laser  using  an  integrated  DFB  laser 
and  electroabsorption  modulator.  The  modulator  facet  was  not  AR  coated  so  the 
reflectivity  from  the  facet  contributes  to  the  phase  change.  The  bandwidth  was  about 
8  GHz  and  the  modulation  efficiency  approximately  10  GHz/mA.  The  linewidth  was 
60  MHz. 

A  number  of  the  Fabry-Perot,  DBR,  and  DFB  FM  lasers  that  have  been  reported 
have  very  high  modulation  efficiency.  This  is  very  desirable  for  microwave  photonic 
links  to  ensure  low  RF  signal  loss  or  even  link  gain.  However,  the  devices  with 
the  highest  efficiency  also  had  some  of  the  worst  linewidths.  None  of  the  linewidths 
reported  were  less  than  1  MHz.  Additional  work  still  needs  to  be  done  to  optimize  FM 
semiconductor  lasers  for  both  low  linewidth  and  high  modulation  efficiency  to  produce 
links  with  low  noise  figures  and  high  dynamic  range.  In  addition,  theoretical  and 
experimental  work  needs  to  be  performed  to  fold  out  the  limitations  of  the  modulation 
nonlinearities  of  multi-section  FM  lasers. 


5.2  Laser  characterization 

The  FM  lasers  we  had  available  for  our  experiments  were  three-section  (gain,  phase, 
grating)  and  four-section  (gain,  phase,  grating,  phase)  distributed  Bragg  reflector 
lasers  designed  by  Bell  Laboratories  and  fabricated  by  Multiplex  Inc.  The  devices  were 
optimized  for  high  modulation  efficiency,  low  linewidth,  and  low  residual  intensity 
modulation. 

I  performed  a  self-heterodyne  measurement  to  characterize  the  linewidth  of  two  of 
the  FM  lasers.  This  technique  was  developed  by  [128],  and  consists  of  heterodyning 
the  laser  with  a  time-delayed  and  AOM  frequency  shifted  version  of  itself.  For  a 
Lorentzian  lineshape,  the  beat  tone  has  a  linewidth  that  is  twice  the  laser’s  linewidth. 
Fig.  5.1  shows  the  experimental  setup  of  the  self- heterodyne  measurement.  Fig. 
5.2  shows  the  80  MHz  beat-tones  from  the  experiment.  Each  laser’s  gain  section  was 
biased  at  200  mA  at  a  temperature  of  20  C.  The  three-section  laser  had  approximately 
a  174  kHz  linewidth,  and  the  four-section  laser  had  a  161  kHz  linewidth. 

I  performed  a  DC  tuning  measurement  to  choose  optimal  bias  points  for  mod¬ 
ulation.  The  lasers  were  gain  biased  at  200  mA,  and  the  phase  section  bias  was 
adjusted.  The  wavelength  was  viewed  on  a  high-resolution  optical  spectrum  analyzer. 
Fig.  5.3  shows  the  tuning  curves  for  the  two  devices.  At  1.3  mA  phase  section  bias, 
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Figure  5.2:  Self  heterodyne  laser  spectrum  measurements  with  Lorentzian  fits. 

3-section  DBR  laser  4-section  DBR  laser 


the  DC  tuning  slope  for  the  3-section  laser  was  22  GHz/mA,  and  the  tuning  slope 
for  the  4-section  laser  was  9  GHz/mA.  At  0.5  mA  phase  section  bias,  the  DC  tuning 
slope  for  the  3-section  laser  was  45  GHz/mA.  The  3-section  device  had  a  much  higher 
efficiency. 

RF  modulation  efficiency  measurements  were  performed  with  a  vector  network 
analyzer.  Fig.  5.4  shows  the  experimental  setup.  A  50  GHz  FSR  (20  pS  delay)  Mach- 
Zehnder  interferometer  was  used  to  discriminate  the  FM.  The  modulation  efficiency 
was  calculated  using  the  gain  expression  (5.1)  corrected  for  the  passive  impedance 
matching  in  the  experiment  and  corrected  for  reflected  modulation  power.  The  mod¬ 
ulation  efficiency  for  both  devices  at  10  MHz  was  about  10%  less  than  at  DC.  At  1.3 
mA  phase  section  bias,  the  bandwidth  of  the  three-section  laser  was  70  MHz,  and 
the  bandwidth  of  the  four-section  laser  was  75  MHz.  At  0.5  mA  phase  section  bias, 
the  bandwidth  of  the  three-section  laser  was  60  MHz.  The  effective  FM  modulation 
efficiency  of  a  lithium  niobate  phase  modulator  is  also  shown  on  the  plot.  For  fre¬ 
quencies  less  than  5  GHz,  the  link  gain  provided  by  the  FM  lasers  is  more  than  an 
order  of  magnitude  better  than  that  of  the  lithium  niobate  phase  modulator. 

Fig.  5.6  shows  the  expected  phase-noise  limited  noise  figure  versus  frequency.  The 
calculation  assumes  a  50  ohm  input  impedance,  and  smoothing  has  been  applied  to 
the  graph.  At  100  MHz  modulation  frequency,  the  modulation  efficiency  of  the  3- 
section  laser  biased  at  0.5  mA  is  high  enough  to  provide  a  6  dB  noise  figure,  but  the 
efficiency  falls  off  rapidly. 


|*S21 


Finally,  I  characterized  the  residual  intensity  modulation  of  the  three-section  DBR 
laser  biased  at  1.3  mA.  A  link  S21  measurement  was  performed  with  the  MZI  to 
measure  the  FM  response,  and  without  the  MZI,  but  at  the  same  photocurrent,  to 
measure  the  IM  response.  In  both  cases,  only  the  phase-section  of  the  device  was 
modulated.  At  100  MHz  modulation  frequency,  the  IM  was  27  dB  lower  than  the  FM. 
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Figure  5.3:  DC  tuning  measurement  of  DBR  laser  phase  sections. 
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Figure  5.4:  FM  modulation  efficiency  experimental  setup. 
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Figure  5.5:  DBR  FM  modulation  efficiency  versus  frequency. 
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Figure  5.6:  Phase-noise  limited  noise  figure  for  FM  DBR  lasers  from  measured  mod¬ 
ulation  efficiency  and  linewidth. 
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Figure  5.7:  Residual  intensity  modulation  measurement  of  DBR  FM  lasers. 
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The  residual  IM  was  sufficiently  small  such  that  it  would  not  degrade  the  linearity 
performance  of  the  link. 

Overall,  the  current-injected,  3-section  DBR  devices  had  high  modulation  effi¬ 
ciency,  low  linewidth,  and  low  residual  IM.  However,  they  had  very  limited  band¬ 
width.  This  is  an  intrinsic  issue  with  current-injection,  as  the  carrier  lifetime  is  very 
long  in  passive  phase-modulation  sections.  Additional  work  needs  to  be  performed  to 
develop  voltage-modulated  devices  in  order  to  increase  the  bandwidth  of  operation  of 
the  FM-DD  microwave  photonic  links. 


5.3  Frequency-modulated  link  with  HR  filter 

The  FM  laser  used  in  the  link-experiments  is  similar  to  the  devices  described  in  the 
previous  section.  It  is  a  three-section  DBR  laser  with  a  45  kHz  linewidth  and  7.5 
GHz  /  mA  modulation  efficiency  at  100  MHz  by  modulating  the  phase  section  of 
the  device.  The  device  was  used  in  conjunction  with  the  HR  discriminator  to  form 
an  FM-DD  link.  Figure  5.8  on  page  73  shows  the  gain  of  the  link  compared  to 
the  PM-DD  link.  For  a  fairly  small  photocurrent  of  5.2  mA  per  detector  (10.5  mA 
total  photocurrent),  we  are  able  to  achieve  positive  link  gain  for  the  FM-DD  link  for 
frequencies  up  to  500  MHz.  Although  the  modulation  efficiency  falls  off  fairly  fast 
because  the  device  is  operated  by  current  injection,  the  link  gain  exceeds  the  lithium 
niobate  phase  modulator  up  past  3  GHz. 

We  measured  the  distortion  of  the  link  and  compare  it  to  the  PM-DD  and  IMDD. 
Figure  5.9  on  page  73  compares  OIP2  and  OIP3  of  the  FM  link  with  PM  link  data 
shown  previously  in  Figure  4.14  on  page  63.  At  its  most  linear  frequency  the  FM 
laser  exceeds  the  OIP3  of  the  Mach  Zehnder  by  5.8  dB,  and  the  link  has  low  second 
order  distortion.  The  link  gain  and  OIP3  degrade  for  higher  modulation  frequencies, 
but  have  significant  improvements  over  IMDD  and  PM-DD  links  at  low  frequencies 
for  both  linearity  and  gain. 


72 


Figure  5.8:  Link  gain  versus  modulation  frequency  for  the  FM  link  versus  the  PM+IIR 
link. 


Figure  5.9:  Distortion  versus  modulation  frequency,  compared  to  the  results  of  the 
PM+IIR  link. 
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5.4  Summary 


FM  lasers  are  a  promising  transmitter  for  microwave  photonic  systems,  mainly  be¬ 
cause  of  their  very  high  modulation  efficiency,  which  exceeds  that  of  the  best  lithium 
niobate  modulators  by  an  order  of  magnitude.  Previous  work  had  not  focused  specif¬ 
ically  on  optimizing  both  laser  linewidth  and  modulation  efficiency  in  order  to  ensure 
low  phase-noise  limited  noise  figures  when  the  links  are  operated  with  high  optical 
power.  We  have  demonstrated  links  that,  within  the  laser’s  modulation  bandwidth, 
provide  noise  figures  <  10  dB,  and  exceed  the  linearity  of  a  MZM  IM-DD  link  when 
using  complementary,  linear-field  PLC  filters  for  demodulation.  With  further  opti¬ 
mization  of  the  laser  modulation  bandwidth  by  using  field-effect  modulation  rather 
than  current-injection,  the  FM-DD  link  architecture  should  be  able  to  provide  mi¬ 
crowave  frequency  links  with  high  dynamic  range,  low  noise  figures,  and  signal  gain. 
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Chapter  6 

Conclusions  and  future  work 


The  work  presented  in  this  dissertation  illustrates  the  benefits  of  phase-modulated  and 
frequency-modulated  microwave  photonic  links  with  direct  detection  (PM/FM-DD) 
as  approaches  to  achieving  microwave  photonic  systems  with  low  noise-figure,  high 
linearity,  and  high  dynamic  range.  By  fully  modeling  the  nonlinear  signal  transmission 
in  the  links,  we  are  able  to  accurately  predict  the  performance  of  a  link  using  particular 
optical  filters  in  the  demodulator.  This  dissertation  has  developed  the  tools  necessary 
for  validating  the  distortion  produced  by  a  particular  filter  designed  to  discriminate 
phase  or  frequency  modulation.  The  close  match  between  experimental  results  and 
simulation,  as  shown  in  Figure  4.10  on  page  60,  was  very  promising.  The  general 
large-signal  and  small-signal  analytical  expressions  for  filtered  PM  and  FM  links,  as 
well  as  our  numerical  model,  should  be  useful  for  future  link  designers. 

On  embarking  on  this  program  of  research,  it  was  to  the  surprise  of  this  author  that 
previous  work  had  relied  on  heuristics  and  inaccurate  Taylor-series  models.  Although 
PM/FM-DD  has  been  actively  explored  since  the  1960’s,  as  discussed  in  Section  2.3, 
we  were  the  first  to  present  complimentary  linear-field  demodulation  as  a  simple 
ideal  model  for  which  to  design  discriminator  filters.  Although  the  early  work  of 
[33]  had  suggested  a  linear-field  transmission  ramp,  they  did  not  identify  the  use  of 
balanced  detection  to  cancel  second  harmonics  or  provide  any  guidance  about  the 
phase  linearity  of  the  filter.  Later  authors  tried  to  linearize  the  slope  of  the  power 
transmission  instead,  which  we  show  is  misguided  (Section  2.6).  It  is  our  hope  that 
the  models  presented  in  this  dissertation  clear  up  remaining  confusion. 

Our  experimental  results  showed  appreciable  improvements  in  the  link  distortion 
over  Mach-Zehnder  modulator  IM-DD  links.  Intermodulation  OIP3  was  improved  by 
over  6  dB  over  the  MZM,  and  improvement  was  shown  for  frequencies  up  through  10 
GHz.  This  research  is  a  first  step,  and  we  believe  we  have  not  reached  the  full  potential 
of  this  approach.  Due  to  difficulties  in  controlling  filter  coefficients,  we  were  not  able 
to  directly  transfer  designs  from  simulation  to  implementation.  However,  the  reverse 
process,  measuring  the  transfer  function  of  a  filter  implementation  and  inserting  it 
into  simulation,  provided  accurate  predictions  of  link  performance.  I  would  like  to  do 
additional  analytical  work  on  HR  flter  design,  and  explore  the  scaling  of  achievable 
SFDR  with  HR  filter  order.  Fig  3.5.  showed  that  we  can  achieve  a  2  dB  increase  in 
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SFDR  per  filter  order  scaling  benefit  from  the  FIR  filter  designs. 

Dne  to  their  diversity  of  applications,  there  is  comparably  more  work  being  per¬ 
formed  in  the  held  of  microwave  photonic  filters  than  in  linearization  of  microwave 
photonic  links.  The  PM/FM-DD  approach  will  benefit  from  synergies  with  that  work. 
Another  photonic-only  approach  to  linear  microwave  photonic  links,  linearized  inten¬ 
sity  modulators,  is  not  as  general,  since  applicable  fabrication  processes  arc  limited  to 
those  with  optically  active  materials,  such  as  lithium  niobate.  In  PM-DD  links,  the 
modulation  is  simple  and  decoupled  from  the  demodulation  process.  Passive  filters 
may  be  easier  to  design  and  optimize  as  a  separate  component. 

Frequency  modulation  of  semiconductor  lasers  is  another  promising  area  of  re¬ 
search.  Links  using  FM  directly  benefit  from  PM-DD  demodulator  designs,  with  the 
added  advantage  of  signal  gain.  There  are  a  number  of  research  topics  to  explore. 
Unanswered  questions  include  theoretical  models  for  and  experimental  characteriza¬ 
tion  of  the  nonlinearity  of  the  frequency  modulation,  and  discovering  the  fundamen¬ 
tal  tradeoffs  between  noise,  modulation  efficiency,  residual  intensity  modulation,  and 
bandwidth  for  the  semiconductor  lasers.  From  my  review  of  the  literature,  it  is  clear 
that  this  design  space  has  not  been  fully  explored.  Chapter  5  shows  that  even  with 
these  unknowns,  we  were  able  to  demonstrate  a  link  with  signal  gain  for  a  small  pho¬ 
tocurrent,  gain  exceeding  a  PM  link  up  through  3  GHz  modulation  frequency,  and 
linearity  exceeding  that  of  an  IM-DD  MZM  link,  if  only  for  a  limited  bandwidth. 

A  future  area  of  research  remains  in  the  area  of  frequency  up-converting  and 
down-converting  links.  Frequency  converting  links,  also  called  IF-to-RF  or  RF-to- 
IF  links,  are  important  in  many  remote  antenna  systems.  Electronic  mixers  produce 
large  spurious  signals,  and  photonic  mixing  is  promising  for  reducing  these  conversion 
spurs.  We  are  looking  at  methods  to  linearize  and  increase  the  dynamic  range  of  these 
frequency-converting  links  by  using  planar- lightwave  circuit  Liters. 
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Appendix  A 
Simulation  code 


In  order  to  design  PM-DD  and  FM-DD  microwave  photonic  systems,  it  is  very  impor¬ 
tant  to  accurately  access  tradeoffs  between  metrics  such  noise,  linearity,  bandwidth, 
power  consumption,  and  cost.  The  following  simulation  code  applies  the  theory  de¬ 
veloped  in  this  dissertation  towards  modeling  these  complex  systems. 


A.l  Small-signal  simulation 

The  small  signal  model  for  a  PM-DD  or  FM-DD  link  is  based  on  the  expressions 
in  Chapter  2.  It  does  not  include  the  effects  of  residual  intensity  modulation,  but 
includes  the  effect  of  the  common-mode  rejection  ratio. 

function  [10,  SIO  ,  Pl_fl  ,  P2_2fl  P3_3fl  ,  P3_2fl_f2]  =  ... 

1  i nk_ss_dis t  or  t  io n  (  fa  ,  fb  ,  Ha,  Hb,  fl  ,  f 2  ,  CMRR,  bl  ,  b2) 
%LINKSS-DISTORTION  Two—  tone  ,  small— signal  model  for  a  PM/FM-DD  link 
%  Inputs  fa,  fb  are  frequency  vectors  of  offsets  from  the  carrier 

%  for  the  filters  ,  with  the  same  length  as  Ha  and  Hb . 

%  Inputs  Ha,  Hb  are  the  complex  filter  transfer  functions 
%  LINK_FFT_DISTORTION  will  linearly  interpolate  Ha  and  Hb 

%  Inputs  fl  ,  f2  are  the  two  modulation  frequencies 

%  Input  CMRR  is  the  common  mode  rejection  ratio,  in  positive  dB 

%  Inputs  bl  ,  b2  are  the  phase  modulation  depths 

% 

%  Output  are  normalized  currents  and  powers 
%  10:  DC  current 

%  SIO:  Summed  DC  current,  used  to  calculate  shot  noise 

%  Pl-fl  •'  Fundamental  at  frequency  fl 

%  P2-2fl  :  Second  harmonic  of  frequency  fl 

%  PSSfl  :  Third  harmonic  of  frequency  fl 

%  P3-2fl-f2 :  IMD3  at  frequency  2fl-f2 

%  Values  assume  1  ohm  load,  1  W  optical  power,  and  1  A/W  efficiency 
%  Multiply  currents  by  R*Popt  to  scale 
%  Multiply  powers  by  Zout*R~2*  Popt  ~2  to  scale 

[X0,  SXO,  Xl_fl  ,  X2_2fl  ,  X3_3fl  ,  X3_2fl_f2]  =  ... 

m  etrics_ss_  distortion  (fa,  fb,  Ha,  Hb,  fl  ,  f  2  ,  CMRR); 

SIO  =  SXO; 
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10  =  X0; 

Pl_fl  =  bl  .'2.*  Xl_fl  ."2/2; 

P2_2f  1  =  bl.~4.*X2_2fl."2/32; 

P3_3fl  =  bl.~6.*X3_3fl. "2/128; 

P3_2fl_f2  =  bl."4.*b2.~2.*X3_2fl_f2.'2/12  8; 

end 

function  [X0,  SXO,  Xl_fl  ,  X2_2fl  ,  X3_3fl  ,  X3_2fl_f2]  =  ... 

metrics_ss_  distortion  (fa,  fb,  Ha,  Hb,  fl  ,  f  2  ,  CMRR) 
%METRICSSS-DISTORTION  Distortion  constants  for  a  PM/FM-DD  link 
%  Inputs  fa,  fb  are  frequency  vectors  of  offsets  from  the  carrier 

%  for  the  filters  ,  with  the  same  length  as  Ha  and  Hb . 

%  Inputs  Ha,  Hb  are  the  complex  filter  transfer  functions 

%  Inputs  fl  ,  f2  are  the  two  modulation  frequencies 

%  Input  CMRR  is  the  common  mode  rejection  ratio  ,  in  positive  dB 

% 

%  Outputs  are  the  distortion  constants 
%  XO:  DC  term 

%  SXO:  Sum  of  XOs ,  used  to  calculate  shot  noise 

%  Xl-fl  :  Fundamental  at  frequency  fl 

%  X2_2fl  :  Second  harmonic  of  frequency  fl 

%  X3-3fl  :  Third  harmonic  of  frequency  fl 

%  X3-2fl_f2 :  IMD3  at  frequency  2fl-f2 

%  DC  component 

XOa  =  interpl(fa, Ha, 0 )  .*conj(interpl(fa  , Ha, 0 )  )  ; 

XOb  =  interpl  ( fb  ,  Hb,  0  )  .  *  conj  ( interpl  ( fb  ,Hb,  0 )  )  ; 

XO  =  func_subCMRR (CMRR, XOa, XOb)  ; 

SXO  =  abs(XOa)  +  abs(XOb)  ; 

%  1st  harmonic 

Xl_fl  =  func_subCMRR (CMRR,  .  .  . 

func_harmonic  ( fa  ,  Ha,  fl  ,  1) 
func_harmonic  (  fb  ,  Hb,  fl  ,  1)); 

%  2nd  harmonic 

X2_2fl  =  func_subCMRR(CMRR,  .  .  . 
func_harmonic  (  fa  ,  Ha,  fl  ,  2) 
func_harmonic  ( fb  ,  Hb,  fl  ,  2)); 

%  3rd  harmonic 

X3_3fl  =  func_subCMRR(CMRR,  . .  . 

func_harmonic  (  fa  ,  Ha,  fl  ,  3)  , .  . . 
func_harmonic  ( fb  ,  Hb,  fl  ,  3)); 

%  IMD3  difference  term 
X3_2fl_f2  =  func_subCMRR (CMRR,  .  .  . 

func_difference  ( fa  ,  Ha,  fl  ,  f  2  ,  2,  1)  ,... 
func_difference(fb,  Hb,  fl  ,  f 2  ,  2,  1)); 

end 

function  X  =  func_harmonic  ( f  ,  H,  fl  ,  n) 
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%  FUNC-HARMONIC  Returns  figure  of  merit  for  harmonic  n*fl 
%  Input  f  is  a  frequency  vector  of  offsets  from  the  carrier 
%  Input  H  is  complex  filter  transfer  function 
%  Input  fl  is  the  modulation  frequency 
%  Input  n  is  the  harmonic 

coeff  =  inline (’2*(  — l)~g./2~n./ factorial (n— g) .  / facto rial(g)  ’  ,  ’n’  ,  ’g  ’ ) ; 
prefactor  =  coeff(n,0); 

X  =  0; 
for  g  =  0:n 

X  =  X  +  coeff (n , g) ./ prefactor . *  ... 

interpl ( f ,H,  ( n— g) . * fl ) . *  ... 

conj(interpl(f,H,  —  g  .  *  fl  ) )  ; 

end 

end 

function  X  =  func_difference(f,  H,  fl,  f  2  ,  n,  p) 

%  FUNC-DIFFERENCE  Returns  figure  of  merit  for  n*fl  —  p*f2 
%  Input  f  is  a  frequency  vector  of  offsets  from  the  carrier 
%  Input  H  is  complex  filter  transfer  function 
%  Inputs  fl  and  f2  are  the  modulation  frequencies 
%  Inputs  n  and  p  are  the  order 

coeff  =  inline  (  [  ’ 2*(  —  1) * (p+g+k) . /2 ~ (n+p) . / f ac t or ial (n— g)  ’ 

’ ./ factorial (g) ./ factorial ( p— k )./ factorial ( k ) ’ ]  ... 

,’nVpVgVk’); 
prefactor  =  coef f (n , p , 0 , 0 ) ; 

X  =  0; 
for  g  =  0:n 

for  k  =  0:p 

X  =  X  +  coeff (n , p , g , k) ./ prefact  or . *  ... 

interpl(f,H,(  n— g)  .  *  f  1  +(— p+k)  .  *  f2  )  .  *  ... 

conj  ( interpl  (  f  ,H,  —  g  .  *  f  1+k  .  *  f2  ) )  ; 

end 

end 

end 

function  X  =  func_sum(f  ,  H,  fl  ,  f 2  ,  n,  p) 

%  FUNCSUM  Returns  figure  of  merit  for  n*fl  +  p*f2 
%  Input  f  is  a  frequency  vector  of  offsets  from  the  carrier 
%  Input  H  is  complex  filter  transfer  function 
%  Inputs  fl  and  f2  are  the  modulation  frequencies 
%  Inputs  n  and  p  are  the  order 

coeff  =  inline  (  [ ’ 2*  ( —  1) ' ( g+k) . / 2  * ( n+p) . / f ac t o r i a  1 (n— g )  ’  ,... 

’ ./ factorial (g) ./ factorial ( p— k )./ factorial ( k ) ’ ]  ... 

,’nVpVgVk’); 
prefactor  =  coef f (n , p , 0 , 0 ) ; 

X  =  0; 
for  g  =  0:n 

for  k  =  0:p 

X  =  X  +  coeff (n , p , g , k) ./ prefact  or . *  ... 

interpl(f,H,( n— g) . *  f 1  +(p— k) . *  f2 ) . *  ... 

conj(interpl(f  ,  H,  —  g  .  *  f  1  — k  .  *  f2  ) )  ; 
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end 


end 

end 

function  X  =  func_subCMRR  (CMRR, Xa,Xb) 

%  FUNCSUBCMRR  Subtracts  two  values  with  a  dB  CMRR 
X  =  abs  (  Xa  -  Xb  +  (Xa  +  Xb)  *10R-CMRR./10)  )  ; 

end 


A. 2  Large-signal  simulation 

The  large  signal  model  uses  the  expansion  from  (2.9). 

function  [10,  SIO  ,  Pl_fl  ,  P2_2fl  ,  P3_3fl  ,  P3_2fl_f2]  =  ... 

1  i nk_l s_d i s t  o r t  i o n  (  fa  ,  fb  ,  Ha,  Hb,  fl  ,  f 2  ,  CMRR  bl  ,  b2) 
%LINK_LS-DISTORTION  Two—  tone  ,  large  signal  model  for  a  PM/FM-DD  link 
%  Inputs  fa,  fb  are  frequency  vectors  of  offsets  from  the  carrier 

%  for  the  filters  ,  with  the  same  length  as  Ha  and  Hb . 

%  Inputs  Ha,  Hb  are  the  complex  filter  transfer  functions 
%  LINK-LS-DISTORTION  will  linearly  interpolate  Ha  and  Hb 

%  Inputs  fl  ,  f2  are  the  two  modulation  frequencies 

%  Input  CMRR  is  the  common  mode  rejection  ratio,  in  positive  dB 

%  Inputs  bl  ,  b2  are  the  phase  modulation  depths 

% 

%  Output  are  normalized  currents  and  powers 
%  10:  DC  current 

%  SIO:  Summed  DC  current,  used  to  calculate  shot  noise 

%  Pl-fl  ■'  Fundamental  at  frequency  fl 

%  P2-2fl  :  Second  harmonic  of  frequency  fl 

%  P3-3fl  :  Third  harmonic  of  frequency  fl 

%  P3-2fl_f2 :  IMD3  at  frequency  2fl-f2 

%  Values  assume  1  ohm  load,  1  W  optical  power,  and  1  A/W  efficiency 
%  Multiply  currents  by  R*Popt  to  scale 
%  Multiply  powers  by  Zout*R~2*  Popt~2  to  scale 

%  DC  component 

XOa  =interpl  (  fa  ,  Ha,  0  )  .*conj(interpl(fa  ,  Ha,  0 )  )  ; 

XOb  =  interpl  ( fb  ,  Hb,  0  )  .  *  conj  ( interpl  ( fb  ,Hb,  0 )  )  ; 

10  =  func  subCMRR (CMRR XOa, XOb)  ; 

SIO  =  abs(XOa)  +  abs(XOb); 

%  1st  harmonic 

Pl_fl  =  func_subCMRR (CMRR,  .  .  . 

func_ls_curr  ( fa  ,  Ha,  fl  ,  f2  ,  bl  ,  b2  ,  1,  0)  ,... 

func_ls_curr(fb,  Hb,  fl  ,  f2  ,  bl  ,  b2 ,  1,  0)).~2/2; 

%  2nd  harmonic 

P2_2f  1  =  func_subCMRR(CMRR,  .  .  . 

func_ls_curr  ( fa  ,  Ha,  fl  ,  f  2  ,  bl  ,  b2  ,  2,  0)  ,... 

func_ls_curr(fb,  Hb,  fl  ,  f 2  ,  bl  ,  b2 ,  2,  0)).'2/2; 

%  3rd  harmonic 
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P3_3fl  =  func_subCMRR(CMRR,  .  .  . 

func_ls_curr  ( fa  ,  Ha,  fl  ,  f 2  ,  bl  ,  b2  , 
func_ls_curr(fb,  Hb,  fl  ,  f  2  ,  bl  ,  b2  , 


3,  0) 

3,  0 ) ) . ” 2 / 2 : 


%  IMD3  difference  term 
P3_2fl_f2  =  func_subCMRR (CMRR,  .  .  . 

func_ls_curr  (  fa  ,  Ha,  fl  ,  f2  ,  bl  ,  b2  ,  2  ,  —1)  , .  .  . 
func_ls_curr  ( fb  ,  Hb,  fl  ,  f2  ,  bl  ,  b2  ,  2,  — l))+2/2; 


end 


function  I  =  func_ls_cur r  ( f  ,  H,  fl  ,  f2  ,  bl  ,  b2  ,  n,  p) 

%  FUNC-LS-Curr  Returns  current  at  frequency  nfl+pf2 
%  Input  f  is  a  frequency  vector  of  offsets  from  the  carrier 
%  Input  H  is  complex  filter  transfer  function 
fmax=  max(  f  )  ; 
fmin=  min(  f  )  ; 

sbs  =  10;  %  Maximum  number  of  sidebands  to  retain 
I  =  0; 

for  g  =  —sbs : sbs+1 

for  k  =  —sbs  :  sbs+1 

valid  =  ( ( g+n)  *  f  1 +(k+p)  *  f2  )  >=fmin  &  ... 

( ( g+n)  *  f  1 +(k+p)  *  f2  )  <=fmax  &  ... 

(g*fl+k*f2)  >=fmin  &  ... 

(g*fl+k*f2)  <=fmax ; 

if  valid 

1  =  1  +  2*  interpl  (  f  ,H,  (  g+n)  *  f  1 +(k+p)  *  f2  )  ... 

.  *  b e s s e  1  j  ( g+n  ,  bl )  .*  besselj  ( k+p  ,  b2 )  ... 

.*conj(interpl(f  ,H,g*  f  1+k*  f2  ) )  ... 

.  *  conj  (besselj  (g  ,  bl )  .  *  besselj  (k  ,  b2 ) )  ; 

end 

end 

end 

end 


A. 3  Numerical  simulation 

The  numerical  model  for  an  FM  link  is  based  on  the  model  in  Figure  3.7  on  page  49. 

function  [10,  SIO  ,  Pl_fl  ,  P2_2fl  ,  P3_3fl  ,  P3_2fl_f2]  =  ... 

link_fft_  distortion  (fa,  fb,  Ha,  Hb,  fl  ,  f  2  ,  CMRR,  bl  ,  b2) 
%LINK  FFT  DISTORTION  Two— tone  ,  numerical  model  for  a  PM/FM-DD  link 
%  Inputs  fa,  fb  are  frequency  vectors  of  offsets  from  the  carrier 

%  for  the  filters  ,  with  the  same  length  as  Ha  and  Hb . 

%  Inputs  Ha,  Hb  are  the  complex  filter  transfer  functions 
%  LINK_FFT_DISTORTION  will  linearly  interpolate  Ha  and  Hb 

%  Inputs  fl  ,  f2  are  the  two  modulation  frequencies 

%  Input  CMRR  is  the  common  mode  rejection  ratio,  in  positive  dB 

%  Inputs  bl  ,  b2  are  the  phase  modulation  depths 

% 

%  Output  are  normalized  currents  and  powers 
%  10:  DC  current 
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%  SIO :  Summed  DC  current,  used  to  calculate  shot  noise 

%  Pl-fl  •'  Fundamental  at  frequency  fl 

%  P2-2fl  :  Second  harmonic  of  frequency  fl 

%  P3-3fl  :  Third  harmonic  of  frequency  fl 

%  P3-2fl-f2 :  IMD3  at  frequency  2fl-f2 

%  Values  assume  1  ohm  load,  1  W  optical  power,  and  1  A/W  efficiency 
%  Multiply  currents  by  R*Popt  to  scale 
%  Multiply  powers  by  Zout*R~2*  Popt'2  to  scale 

%%  FFT  setup  %% 

%  Make  sure  discretation  captures  modulation  frequencies 
fstep  =  min  ( [  abs  (  f2  —  f  1  )  ,  fl,  f2  ]  )  ; 

flmod  =  round  (fl  ./  fstep)  .*  fstep  ;  f2mod  =  round  (  f  2  .  /  fstep)  .*  fstep  ; 

%  Maximum  sampling  frequency  is  limited  by  the  spectral  domain  provided 
fmax  =  min  ( [  abs  (max(  fa  ) )  ,  abs  ( min(  fa  ) )  ,  abs  (max(  fb  ) )  ,  abs  (min  (  fb  ) )  ] )  ; 

%  Number  of  points  (N  is  odd  to  include  DC) 

N  =  2*  floor  ( fmax ./ fstep  ) +1; 

%  Sampling  frequency 
fS  =  fstep*N; 

%  One-sided  frequency  vector 
FIs  =  (0:N/2-l/2)*fS/N; 

%  Two-sided  frequency  vector 
F2s  =  (— N/2  +  l/2:N/2  — l/2)*fS/N; 

%  Time  vector 
t  =  ( 0  :N— 1)  /  fS  ; 

%  Resample  the  filters 
Ha2s  =  interpl  (  fa  ,Ha,  F2s)  ; 

Hb2s  =  interpl  (  fb  ,Hb,  F2s)  ; 

%%  Modulation  %% 

%  Instantaneous  phase  deviation  signal 

XinT  =  bl *cos  (2*  pi *flmod*  t )  +  b2*cos  (2*  pi *f2mod*  t )  ; 

%  Modulated  optical  signal  centered  at  fc  ,  1W  of  power 
EmodT  =  sqrt  (2)  ,*exp(l  i*XinT)  ; 


%%  Filtering  %% 

EmodF  =  fftshift  (  fft  (EmodT, N)  )/N;  %  Convert  to  the 
EoutaF  =  EmodF.*  Ha2s;  EoutbF  =  EmodF.  *  Hb2s  ;  %  Split 
EoutaT  =  if  ft  (  i  f  f  t  s  h  i  f  t  (EoutaF)  ,N)  *N;  %  Convert  to 

EoutbT  =  ifft  (  ifftshift  (EoutbF)  ,N)*N;  %  Convert  to 


frequency  domain 
and  Filter 
time  domain 
time  domain 


%%  Detection  %% 

%  Convert  back  to  electrical  domain 
IoutaT  =  EoutaT  .*  conj  (EoutaT) /2; 

IoutbT  =  EoutbT  .*  conj  (EoutbT)  / 2 ; 

IoutT  =  IoutaT  —  IoutbT  +  (IoutaT  +  IoutbT  )  *10' (-CMRR.  /  1 0 )  ; 
10  =  mean(  IoutT); 

SIO  =  abs  (mean(  IoutaT )  )+abs  (mean(  IoutbT  ))  ; 


%  This  gives  the  RMS  power:  multiply  by  2  for  two  sided 
%  spectrum ,  square  it  ,  then  divide  by  2  for  RMS 
%  The  DC  power  should  not  be  multiplied  by  2. 

IoutF  =  fftshift  (  fft  (IoutT  ,N)  )/N; 

PoutF  =  2*abs  ( IoutF  (N/2  +  l/2:end) )  . '  2  ; 
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PoutF(l)  =  PoutF(l)  / 2 ; 


%%  Analysis  %% 

Pl_fl  =  interpl  (FIs  ,  PoutF  ,  flmod)  ; 

P2_2fl  =  interpl  (FIs  ,  PoutF  ,2*  flmod)  ; 

P3_3fl  =  interpl  (FIs  ,  PoutF  ,3*  flmod)  ; 

P3_2fl_f2  =  interpl  ( FIs  ,  PoutF  , 2 *  flmod— f2mod )  ; 

end 


A. 4  Link  response 

This  program  will  output  the  response  of  the  link,  for  a  given  input  power,  for  any  of 
the  three  models. 

function  [Idc,  Side,  Pl_fl  ,  P2_2fl  ,  P3_3fl  ,  P3_2fl_f2]  =  ... 

link_response  ( fa  ,  fb  ,  Ha,  Hb,  fl  ,  f  2  ,  CMRR,  PI,  P2 ,  param) 
%LINK_RESPONSE  Two-tone  model  for  a  PM/FM-DD  link 
%  Inputs  fa,  fb  are  frequency  vectors  of  offsets  from  the  carrier 
%  for  the  filters  ,  with  the  same  length  as  Ha  and  Hb . 

%  Inputs  Ha,  Hb  are  the  complex  filter  transfer  functions 

%  Inputs  fl  ,  f2  are  the  two  modulation  frequencies 

%  Input  CMRR  is  the  common  mode  rejection  ratio,  in  positive  dB 

%  Inputs  PI,  P2  are  the  modulation  power  at  each  frequency 

%  Input  param  is  a  cell  array  containing 
%  modtype:  ’ PM ’  or  ’FM’ 

%  model:  ’  smallsignal  ’ ,  ’  largesignal  or  ’  fft 

%  Vpi:  modulator  halfwave  voltage  for  PM 

%  eta:  modulation  efficiency  Hz/A  for  FM 

%  Side:  desired  total  detector  power  to  scale  to 

%  Zin :  input  impedance  magnitude 

%  Zout:  output  impedance  magnitude 

%  Iw :  linewidth  of  laser 

% 

%  Output  are  currents  and  powers 
%  Id:  DC  current 

%  SIO :  Summed  DC  current,  used  to  calculate  shot  noise 

%  Pl-fl  •'  Fundamental  at  frequency  fl 

%  P2-2fl  :  Second  harmonic  of  frequency  fl 

%  P3-3fl  :  Third  harmonic  of  frequency  fl 

%  P3-2fl_f2 :  IMD3  at  frequency  2fl-f2 

[  modtype  ,  model  ,  Vpi ,  eta  ,  Side  ,  Zin  ,  Zout  ,  lw  ]  =  param  {:}  ; 

II  =  sqrt  (2*P1/Zin )  ;  12  =  sqrt  (2*P2/Zin)  ; 

VI  =  sqrt (2*P1* Zin ) ;  V2  =  sqrt (2*P1* Zin ) ; 

switch  modtype 
case  ’PM’ 

bl  =  pi*Vl./Vpi; 
b2  =  pi*V2./Vpi; 
case  ’FM’ 

bl  =  eta*Il ./ fl  ; 
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end 


b2  =  eta  *  12 . / f2  ; 


switch  model 

case  ’  smalls ignal  ’ 

[10,  SIO  ,  Pl_fl  ,  P2_2fl  ,  P3_3fl  ,  P3_2fl_f2 

link_ss_distortion(fa,  fb,  Ha,  Hb,  fl  , 
case  ’largesignal  ’ 

[10,  SIO,  Pl_fl  ,  P2_2fl  ,  P3_3fl  ,  P3_2fl_f2 

link_ls_distortion(£a,  fb  ,  Ha,  Hb,  fl  , 
case  ’ fft  ’ 

[10,  SIO,  Pl_fl  ,  P2_2fl  ,  P3_3fl  ,  P3_2fl_f2 

link_fft_  distortion  (fa,  fb,  Ha,  Hb,  fl 

end 


f2  ,  CMRR,  bl ,  b2)  ; 
f2  ,  CMRR,  bl ,  b2)  ; 
f2  ,  CMRR,  bl  ,  b2)  ; 


RPoptScale  =  Side.  /SIO; 

Idc  =  10  .*  RPoptScale  ; 

Pl_fl  =  P 1  _f  1  .*  Zout  .*  RPoptScale  .  ~  2 ; 

P2_2f  1  =  P2_2f  1  .  *  Zout .  *  RPoptScale  .  ~  2  ; 

P3_3fl  =  P3_3fl.*  Zout.*  RPoptScale.  *2; 

P3_2fl_f2  =  P3_2fl_f2.*Zout.*  RPoptScale  . 


2; 


end 


A. 5  Link  metrics 

This  program  will  calculate  the  small  signal  metrics  of  a  given  link. 

function  [  OIP2dBm ,  OIP3dBm ,  GdB ,  NdBm,  NFdB ,  SFDRdB  ]  =  ... 

link_metrics  ( fa  ,  fb  ,  Ha,  Hb,  fl  ,  f 2  ,  CMRR,  param) 

%LINKJMETRICS  Small  signal  metrics  for  a  PM/FM-DD  link 
%  Inputs  fa,  fb  are  frequency  vectors  of  offsets  from  the  carrier 
%  for  the  filters  ,  with  the  same  length  as  Ha  and  Hb . 

%  Inputs  Ha,  Hb  are  the  complex  filter  transfer  functions 

%  Inputs  fl  ,  f2  are  the  two  modulation  frequencies 

%  Input  CMRR  is  the  common  mode  rejection  ratio,  in  positive  dB 

%  Inputs  PI,  P2  are  the  modulation  power  at  each  frequency 

%  Input  param  is  a  cell  array  containing 
%  modtype:  ’PM’  or  ’FM’ 

%  model:  ’  smallsignal  ’ ,  ' largesignal  ’ ,  or  ’fft 

%  Vpi:  modulator  halfwave  voltage  for  PM 

%  eta:  modulation  efficiency  Hz/A  for  FM 

%  Side:  desired  total  detector  power  to  scale  to 

%  Zin :  input  impedance  magnitude 

%  Zout:  output  impedance  magnitude 

%  Iw :  linewidth  of  laser 

%  Output  are 

%  OIP2dBm:  OIP2  in  dBm 

%  OIPSdBm:  OIP3  in  dBm 

%  GdB:  Gain  in  dB 

%  NdBm:  Receiver  noise  in  dBm 

%  NFdB:  Noise  figure  in  dB 
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SFDRdB:  SFDR  in  dB  in  1  Hz  bandwidth 


q  =  1 . 6 e  — 19;  %  Electron  charge  (C) 

I<T  =  4e  — 18;  %  Thermal  noise  at  room  temperature  (W) 

[  modtype  ,  model  ,  Vpi ,  eta  ,  Side  ,  Zin  ,  Zout  ,  lw  ]  =  par  am  {  :  }  ; 

[XO,  SXO,  Xl_fl  ,  X2_2fl  ,  X3_3fl  ,  X3_2fl_f2]  =  ... 

metrics_ss_distortion  (fa  ,  fb  ,  Ha/sqrt(2),  Hb/sqrt(2),  fl  , 
RPopt  =  Side. /SXO; 

Ssn  =  2*q* Side* Zout ; 

Spn  =  Zout  *  RPopt . ”  2  *  lw  .  *  Xl_fl  .  ~  2/ pi .  /  f  1  .  ~  2  ; 

Stn  =  KT; 


switch  modtype 
case  PM’ 

G  =  Zin*  Zout  *  ( Xl_fl  .  *  pi .  *  RPopt/ Vpi )  .  *  2  ; 
case  ’EM’ 

G  =  Zin/Zout  *  ( Xl_fl  .  *  eta  .  *  RPopt .  /  f  1 )  +  2  ; 

end 


OIP2dBm  = 

OIP3dBm  = 

GdB 

NFdB 

NdBm  = 
SFDRdB  = 

end 


10*  loglO  (8*  Zout  *RPopt . '  2 )  +  40*  loglO  ( Xl_fl )  — 
20*logl0(X2_2fl)  +  30; 

10*  loglO  (4*  Zout  *RPopt . '  2 )  +  30*  loglO  ( Xl_fl )  — 
10* loglO  (X3_2fl_f2)  +  30; 

10*  loglO  (G)  ; 

10*loglO(l  +  l./G+(Spn+Ssn)  ./ (G*KT) )  ; 

10*  loglO  (  Spn+Ssn+Stn )  +30; 

2/3*(OIP3dBm— 30— GdB+174— NFdB)  ; 


f2  ,  CMRR) 
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